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ABSTRACT 


This  thesis  deals  with  free  surface  flows  of  an  ideal  fluid  in 
a  gravity  field.  The  work  is  restricted  to  two  dimensions  and  only  irrota- 
tional  flows  are  considered.  Lewy’s  proof  that  the  free  surface  of  such 
flows  is  an  analytic  curve  is  given  in  Chapter  Two,  while  Chapter  Three 
gives  an  account  of  some  methods,  based  on  Lewy’s  result,  for  obtaining 
the  flow  corresponding  to  a  particular  free  surface.  Some  connections 
between  the  various  methods  and  a  re-derivation  of  one  of  them  are  given 
in  Chapter  Four.  Chapter  Five  contains  some  examples  of  flows  of  this 
type;  most  of  these  involve  the  flow  over  a  sharp  corner  of  a  stream 
which  is  uniform  and  of  finite  depth  at  great  distances  upstream.  Such  a 
corner  exists  in  the  flows  corresponding  to  a  wide  class  of  free  surface 
curves,  and  some  discussion  of  this  phenomenon  occurs. 
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CHAPTER  ONE 

INTRODUCTION 


We  shall  concern  ourselves  here  with  the  flow  of  a  fluid  with  a 
free  surface  under  the  action  of  gravity.  We  shall  assume  that  the  fluid 
under  consideration  is  non-viscous  and  incompressible,  and  that  the  motion 
is  both  steady  and  irrotational.  By  steady  flow  we  mean  that  the  situation 
at  a  fixed  point  in  space  does  not  change  as  time  goes  on;  by  irrotational 
flow  we  mean  that  if  the  velocity  at  a  point  in  the  fluid  Is  expressed  as 
a  vector  function  of  position,  then  this  vector  has  zero  curl.  We  shall 
also  assume  that  the  motion  is  two-dimensional;  that  is  to  say  we  shall 
assume  that,  by  a  suitable  choice  of  cartesian  axes,  we  can  make  every 
physical  variable  independent  of  one  coordinate.  It  will  suffice,  then,  to 
discuss  the  motion  in  the  plane  of  the  other  two  coordinates.  The  term 
free  surface  or  free  streamline  will  mean  any  streamline  along  which  the 
pressure  does  not  vary.  This  condition  is  met  on  the  surface  of  a  fluid 
open  to  the  air.  A  gravity  flow  will  be  taken  to  mean  a  flow  where  the. 
effects  of  gravity  are  taken  into  account.  If  these  effects  are  neglected 
we  shall  speak  of  a  gravity-free  flow. 

We  choose  a  coordinate  system,  then,  with  the  y-axis  pointing  up 
so  that  gravity  acts  in  the  negative  y  direction  Denote  by  u(x,y)  and 
v(x,y)  the  x-  and  y-components  of  the  fluid  velocity  at  the  point  (x,y). 
The  conditions  of  incompressibility  and  irrotational ity  imply  that,  inside 
the  fluid,  the  quantities  u(x,y)  and  -v(x,y)  satisfy  the  Cauchy -Riemann 
equations.  Hence  the  complex  velocity  w  of  the  fluid  given  by 

w  =  u  -  iv 
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is  an  analytic  function  of  the  complex  variable 

z  =  x  +  i  y 

inside  the  fluid.  The  speed  of  the  fluid  is  denoted  by  V  and  is  given 
by 

v  -  M  • 

Since  u>  is  analytic,  so  is  its  integral,  and  so  the  complex  potential 
f unc  t ion 

f  =  <p(x,y)  +  i  ^(x,y) 

defined  by 

df 

dz  -  W 

is  also  an  analytic  function  inside  the  fluid.  This  function  is  arbitrary 
up  to  a  complex  constant  and  it  is  not  hard  to  show  that  the  functions 
cp(x,y)  and  ^(x, y)  are  respectively  the  velocity  potential  and  the  stream 
function  for  the  flow.  That  is  to  say  the  gradient  of  cp  is  just  the 
velocity  vector  function,  and  t'he  curves  ^  =  constant  are  the  streamlines 
of  the  flow.  We  choose  the  arbitrary  constant  in  ijj  so  that  the  free 
surface  of  the  fluid  is  given  by  ^  =  0. 

Since  the.  fluid  is  incompressible  and  irrctat::  Mi.al  then  Bernoulli5 s 

equation  gives  that 

\  V2  +  gy  +  ^ 

is  constant;  throughout  the  fluid.  In  this  expression  p  denotes  the 
pressure  while  p  is  the  fluid  density.  The  condition  that  distinguishes 
the  free  surface  is  that  the  pressure  does  not  vary  there.  Hence,  if  U  is 
a  constant  with  the  dimension  of  velocity,  then 


• ;  1  ' 

'  ■.  '  .  U  v 
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V2  +  2gy  =  U2 

is  Bernoulli's  equation  on  i]j  =  0.  By  suitable  choice  of  origin  the  constant 
U  can  be  chosen  equal  to  zero;  this  is  the  situation  in  Chapters  Two  and 
Three.  In  Chapters  Four  and  Five  the  more  general  situation  is  dealt  with. 

From  a  practical  point  of  view  the  most  interesting  problem  is 
that  in  which  the  shape  of  the  bed  and  the  upstream  conditions  are  given, 
and  the  free  streamline  is  to  be  determined.  This  problem  is  difficult 
because  it  is  Impossible  to  compute  V  on  the  bed  using  Bernoulli" s  equa 
tion  since  the  pressure  there  is  not  known.  On  the  other  hand,  the  press¬ 
ure  is  constant  on  the  free  streamline  and  so  the  velocity  is  known  if  the 
curve  is  given.  This  makes  it  somewhat  easier  to  determine  a  flow  consistent 
with  that  free  streamline.  The  methods  dealt  with  here  all  solve  this  latter 
problem;  they  assume  that  the  free  streamline  is  a  known  curve  and  proceed 
from  this  using  Bernoulli's  equation. 

From  our  point  of  view,  the  most  important  property  of  the  flows 
which  we  are  considering  is  that  the  free  streamline  is  an  analytic  curve. 
This  was  first  proved  by  Lewy  [1]  in  1952  and  the  proof  is  presented  in 
Chapter  Two,  Using  this  fact,  analytic  parametric  expressions  can  be  found 
for  f  and  z  which  are  valid  on  the  free  streamline.  These  ax©  assent"* 
lally  re  “•statements  of  Bernoulli’s  equation  and  the  equation  of  the  free 
streamline,  and  we  can  arrange  It  so  that  the  parameter  involved  is  real 
on  m  0.  Since  these  are  analytic  relations,  thair  power-series-  expans¬ 
ions  can  be  used  to  define  analytic  functions  for  complex  values  of  the 
parameter,  and  these  in  turn  determine  a  flow  consistent  with  the  given  free 
streamline.  This  process  has  been  carried  out  by  Lewy  [2]  and  John  [3] 
using  quite  different  formulations  of  the  problem.  Milne- Thomson  [4]  extends 
the  results  of  John.  These  methods  are  presented  in  Chapter  Three  and,  in 
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Chapter  Four,  Lewy ' s  results  are  re-derived  in  a  different  fashion  and  some 
connections  are  made  between  these  results  and  those  of  John  and  Milne- 
Thomson. 


The  flows  constructed  in  this  way  cannot  in  general  be  extended 
away  from  the  free  streamline  past  a  certain  point.  This  point  corresponds 
to  a  singularity  in  the  complex  velocity  and  its  location  depends  on  the 
shape  of  the  free  streamline.  Any  streamline  above  or  through  this  point 
can  of  course  be  chosen  as  a  bed.  The  streamline  through  the  singularity 
often  has  a  sharp  corner  at  that  point?  the  angle  and  orientation  of  this 
corner  depend  on  the  flux  of  fluid  down  the  stream  as  well  as  on  the  shape 
of  the  free  streamline.  This  dependence  is  discussed  in  Chapter  Five,  and 
some  illustrative  examples  are  presented  there.  The  first  five  of  these 
are  uniform  flows  of  finite  depth  far  upstream,  and  they  all  have  a  corner 

2jt 

of  angle  “  .  This  does  not  mean  that  other  angles  are  impossible  and 

3 

Example  6,  while  it  does  not  have  the  same  upstream  behavior,  exhibits  how 

different  angles  can  occur.  In  addition  there  exist  flows  which  have  no 

such  singularity?  the  flow  with  a  straight  horizontal  free  streamline  and 

the  flottf  in  the  example  at  the  end  of  Chapter  Four  are  examples  of  this. 

However  the  latter  flow  has  a  corner  of  angle  -r-  in  the  free  streamline. 
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CHAPTER  TWO 

THE  ANALYTIC ITY  OF  THE  FREE  STREAMLINE 


The  method  used  here  was  given  by  Hans  Lewy  in  1952.  [1].  It  is 
used  to  establish  that  the  free  streamline  in  a  flow  such  as  we  are  consider¬ 
ing  is  an  analytic  curve.  Lewy's  result  is  more  general;  the  relation  which 
he  imposes  on  the  boundary  of  the  region  of  concern  is  an  arbitrary  analytic 
function  whereas  here  we  have  the  explicit  relation  embodied  in  Bernoulli" s 
equation.  The  method,  however,  is  the  same. 

Bernoulli's  equation  on  ip  =  0  reads: 

lff|2  -  -  2gy  .  t  =  0  (2.01) 

Here  f  =  cp  +  tip  =  f  ( z  )  is  an  analytic  function  of  z  =  x  +  iy  in  the 
fluid,  i.e.  in  the  region  if/  <  0.  At  a  point  where  J”j  =  v  /  0,  f  =  f(z) 
can  be  inverted  to  give: 

z  =  F(f)  =  x(q ),ip)  +  iy(q>»^)  =  F(<P>^)  (2.02) 

which  is  an  analytic  function  of  f  =  cp  +  i3p  in  some  neighborhood  of  a 
point  where  V2>  0.  Using  equation  (2.01)  we  wish  to  show  that  F(f)  can 
be  extended  analytically  across  the  line  ip  =  0  wherever  >  0. 

Now  where  2  >  0  we  have  3^  =  [4§]  ^  =  [77]  ^  •  Hence 

dz  Ldf  udf 

equation  (2.01)  becomes: 

1  =  -  2gy  |“|2  .  ip  =  0,  V2  >  0.  (2.03) 

Consider  the  function 


6(f)  =  F(f)  . 


(2.0l*) 


■  v  c'  -  y  ?  1  ■  J 
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On  ip  =  0  f  is  real  and  so  G(f)  =  F(f).  Hence  we  have  that  2iy  =  F  -  F 
=  F  -  G.  Moreover 


,dF|2  _  dF  ,dF\  _  dF  dF  _  dF  dG 
'df1  "  df  df '  ~  df  df  "  df  df 

on  ip  =  0,  and  so  equation  (2.03)  can  be  written 


1 


ig(F  -  G) 


dF 

df 


dG 

df 


It  follows  that 

~  =  i[g  ~  (G-F)]"1  =  H(G,f).  ip  =  0.  (2.05) 

At  points  where  V2  >  0,  F(f)  is  analytic  and,  in  addition,  equation 
(2.01)  ensures  that  y,  and  hence  G-F,  is  non-zero.  Hence  H(G,f)  is 
an  analytic  function  of  both  its  arguments  where  Vs  >  0. 


Inside  the  fluid  where  F  and  G  are  both  analytic  functions 
of  f,  we  can  write: 


If  =  If  G('’’^  ■  • 1  % 


df  =  =  •  1  %  • 


(2.06) 

(2.07) 


In  order  to  establish  the  analyticity  of  the  free  streamline  we 

must  make  some  assumptions  about  x  and  y.  These  are  that  x(cp,^), 

dx 

y(op,^)  and  ^  exist  and  are  continuous  for  ip  <  0.  Equation  (2.02) 
implies  that  x  and  y  are  conjugate  harmonics  with  continuous  derivatives 
inside  the  fluid.  Our  assumptions  together  with  the  mean  value  theorem 
imply  that  ^  exists  and  equals  for  ip  <  0.  (See  Appendix  1).  Also 

equations  (2.O7)  and  (2.02)  together  with  (2.03)  imply: 


and 
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=  -  +  <§f)2] 

(2.08) 

(2.09) 

on  ip  =  0.  Hence,  also, 
ip  <  0. 


exist  and  are  continuous  in  the  region 


Now,  at  the  origin  equation  (2.04)  yields 


G(0)  =  F(0)  .  (2.10) 


This  condition  together  with  the  differential  equation  (2.05)  can  be  shown 
(see  Appendix  2)  to  ensure  the  existence  and  uniqueness  of  an  analytic 
function  G(f)  inside  a  region  where  certain  conditions  on  H(G,f)  are 
satisfied.  These  conditions  are: 


|H(G,f)|  <  M 

and  l||!  <  N  * 

where  M  and  N  are  constants. 


(2.11) 

(2.12) 


To  show  that  equation  (2.1l)  is  verified  we  remember  that,  since 
|ff|  =  iff I  =  !/v  >  equation  (2.05)  gives 

|H(G,£)|  =  1/g  |£|  |G-F |  =  v/g  |G-F|  . 

Moreover  equation  (2.10)  yields  |G(o)  -  F(o) J  =  |F(o)  -  F(o)j  =  2 1 y J . 
Where  V2  >  0  equation  (2.0l)  ensures  |y|  >0  and  so  | G-F j  >0  at 
the  origin.  Finally,  since  y  is  assumed  to  be  continuous  up  to,  and 
including  the  line  ip  =  0,  we  have  that  |G-F|  is  bounded  away  from  zero 
in  some  small  neighborhood  ip  <  0  of  the  origin.  This  means  that  the  con 
dition  (2.11)  is  met. 
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Using  equation  (2.05)  once  more  we  obtain 


dF 

df 


(G  -  F)2  . 


8h 

This  yields  |^|  <  N  in  the  same  region,  and  for  the  same  reasons,  as 
above.  Hence  condition  (2.12)  is  satisfied. 


This  means  that  equation  (2.05)  can  be  solved  for  a  function  G(f) 
satisfying  equation  (2.10)  which  exists,  is  unique,  and  is  analytic  in  a 
small  neighborhood  ip  <  0  of  the  origin.  Its  analyticity  on  ip  =  0  will 
be  shown  below.  It  is  worthwhile  to  note  that  by  appropriate  choice  of  the 
arbitrary  constant  in  the  velocity  potential  cp  the  origin  can  be  taken  at 
any  point  on  the  free  streamline,  and  so  the  above  argument  applies  at  any 
point  where  V2>  0. 


We  now  introduce  a  function  G(q>,ip)  by  the  following  construction: 

6(0,0)  =  F(0,0)  ,  (2.1?) 

^  0(0 ,*)  -i[g  ^  (G(O.l(')-F(O^))]'1  :  H^G.V-)  (2.1!*) 

on  q)  =  0,  ip  <  0,  and 

6(q>,^0)  =  l[g  ^  F(cp.^0)  (G(q>.^0)-F(cp,^0)))_1  S  H^G.Cp,^) 

(2.15) 

on  ip  =  ip  <  0.  Equation  (2.14)  can  be  considered  as  an  ordinary  differ¬ 
ential  equation  for  G(0,ip).  Together  with  the  condition  (2.13)  it  deter¬ 
mines  a  function  G(0 ,ip)  which  exists,  is  unique  and  continuous  (see 
Appendix  2)  in  a  region  ip  <  0  of  the  negative  ip  axis  where  conditions 
comparable  to  (2.11)  and  (2.12)  are  met  by  H^(G,^).  The  proof  that  these 
conditions  are  satisfied  parallels  that  given  for  H(G,f)  above.  The 
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continuity  follows  since  F 
by  hypothesis. 


and 


6f 

cty/ 


are  continuous  up  to  the  free  surface 


Having  thus  constructed  G(0 yp) ,  we  use  this  boundary  value  with 
equation  (2.15).  Pq  must,  of  course,  be  chosen  in  the  region  of  existence 
of  G(0  ,if/) .  Again  (see  Appendix  2)  the  function  G(cp  ,p^)  defined  in  this 
way  exists  and  is  unique  in  a  region  where  conditions  analogous  to  (2.1l) 

and  (2.12)  are  satisfied  by  K^GjCp >P q)  •  The  region  here  is  p  =  p^  <  0, 

I<P|  <  e  for  some  e  .  The  function  G(0,p^)  is  continuous  for  f0<° 
as  is  Hg(G,cp,^ q)  by  hypothesis.  Hence  (see  Appendix  2)  the  solution 
G(cP,^q)  is  continuous  in  both  cp  and  the  parameter  p^  for  p  <  0.  It 
follows  that  the  function  G(q>,p)  tends  uniformly  to  G(cp,o)  as  p  -+0. 

We  have  that  G(cp,0)  is  a  solution  of  the  differential  equation 
(2.15)  with  initial  value  G(0,0)  =  F(0,0)  .  However  a  solution  satisfying 
this  condition  is  certainly  G(cp,0)  =  F(cp,0),  since  this  substitution 

changes  (2.15)  into  (2.09).  Hence  uniqueness  requires 

G(q>,0)  =  F(cp,0)  .  (2.16) 

Now  if  we  specialize  derivatives  of  G(f)  as  constructed  from 
equations  (2.05)  and  (2.10)  by  the  relations  (2.06),  then  on  cp  =  0,  (2.05) 
becomes  (2.14),  and  on  p  =  p  <  0,  (2.05)  becomes  (2.15).  The  uniqueness 

of  all  these  solutions  then  requires  that  G(f)  =  G(cp  ,p) .  Moreover,  this 
shows  that  the  analytic  function  G(f)  remains  continuous  up  to  the  free 
streamline  as  p  ->-0.  Equation  (2.16)  now  becomes 


G(cp  +  10)  =  F(cp  +  i0)  . 


(2.17) 


■ 


{(.  + 
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We  now  define  for  f  >  0 

F(q>  +  if)  =  G(cp  -  if).  (2.18) 

This  is  clearly  analytic  for  f  >  0  since  G(cp  +  if)  is  analytic  for 
f  <  0.  This  definition  of  the  analytic  function  F(f)  for  f  >  0  is  the 
same  as  our  previous  one  on  the  free  streamline  by  equation  (2.I7),  and  so 
provides  the  required  analytic  extension  of  F(f)  across  the  f- axis. 

This  means  that  F(f)  is  analytic  on  the  free  streamline  and  so  x(cp,^/) 
and  y(cp  ,f)  are. 

Since  y  is  analytic  on  f  =  0,  equation  (2.01)  shows  that  7 
is.  Moreover,  if  we  denote  distance  along  the  free  streamline  by  s,  then 
this  is  analytic  since 

S  =  I  ^  +  dX  ' 

ds 

Finally,  since  -r-  =  V  where  t  is  time,  we  have 


and  so  this  quantity  too  is  analytic.  The  fact  that  all  these  quantities 
are  analytic  functions  of  each  other  will  be  used  extensively  in  the 
following. 
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CHAPTER  THREE 

METHODS  OF  OBTAINING  FLOWS 


In  this  chapter  we  outline  three  methods  of  obtaining  flows  when 
the  equation  of  the  free  streamline  is  known.  The  first  two  are  independent 
approaches  to  the  problem  while  the  third  is  an  extension  of  the  second. 


Section  5.1  The  Method  of  Lewy. 


The  following  method  appeared  in  1952  in  a  paper  by  Hans  Lewy  [2], 
It  makes  use  of  the  fact  that  the  free  streamline  is  an  analytic  curve  to 
develop  a  way  of  obtaining  a  flow  having  a  given  free  streamline.  The  nota¬ 
tion  here  follows  that  of  Chapter  Two. 


Whenever  |~j2  =  V2  ^  0,  the  equation  f  =  f(z)  has  an  inverse 
and  so,  as  in  equations  (2.02)  and  (2.0^),  we  can  define  functions  F(f)  and 
G(f)  as 


z  =  z(f)  =  F(f) 


(3.01) 


and  _ 

G(£)  =  F(f) 


(3.02) 


On  the  free  streamline,  if/  =  0  and  so 

y  =  -  |[F(£)  -  G(£)] 

a  i  •  I dF | p  dF  dG  .  _ 

Also,  since  l^f I  =  df  df  °n  ^  =  °» 

(2.03)  as 


f  is  real.  This  gives 

f  =  0  (3.03) 

we  can  write  the  Bernoulli  equation 


1  =  lg(F  -  G)  =  0  (3.04) 


We  now  define  a  quantity 

*(£)  =  |[F(£)  -  G(f )] 


(3-05) 
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This  is  a  real  quantity  on  if/  =  0  and,  comparing  equation  (3.O3)  and  (3„05)} 
we  see  that  on  the  free  streamline 


A  =  "  y  .  ^  =  0  {3.06) 

Using  equation  (3-05)  in  equation  (3.04)  yields 

1  •  2gA  H  df  •  *  ‘  0  (5-07) 

We  now  solve  equations  (3. 05)  and  (3-07)  for  F(f)  in  terms  of  the  function 

A(f).  If  we  differentiate  equation  (3.O5)  with  respect  to  f,  solve  the 
dG 

result  for  —  and  substitute  this  into  equation  (3-07)  we  get,  after  some 
simplification  that 


This  gives 


*2 


+  2i 


dA  dF 
df  df 


1 

2gA 


0 


dF  dA  r  1 

df  ”  df  ™  *-2gA 


Mf;  J 


On  the  free  streamline  A  and  f  are  both  real.  Hence,  since  F  =  x+iy* 
the  real  part  of  this  is 


dx  _  ,  1 

dcp  ~  -  ®-2gA 


2J2 


provided  >  (*)2  . 

positive  root  and  so 


In  order  that  x  increase  with  cp 


F(f)  =  z(f)  =  --  iA  + 


(||)2]2  df 


we  choose  the 


(3-08) 


When  A(f)  is  known,  equation  (3.O8)  gives  z  as  a  function  of  f  =  cp+i^t 
In  theory  this  can  be  inverted  where  V2  is  not  singular  to  yield  f  as  a 
function  of  z  and  hence  give  the  velocity  potential  and  stream  function,, 


If  we  remember  that  ”2  =  |~^|2  =  df  df  ^  then  equation  (3-07) 

becomes 
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V2  =  2gA  . 


ip  =  0 


(3-09) 


This  gives  some  conditions  on  the  function  A(f).  Equation  (3.06)  implies 
that  A(f)  is  real  on  ip  --  0,  that  is  it  is  real  for  real  arguments.  In 
addition  equation  (3.O9)  gives 


and 


A  = 

0 

if 

< 

i\j 

11 

0  , 

A  > 

0 

if 

A 

> 

0  . 

ip  =  0  (3.10) 


Now  z(f)  =  x(f)  +  iy(f)  where  x(f)  and  y(f)  are  real  for  real  arguments. 

tt  dz  dx  .  dy 

Hence  -77  =  —?  +  1  —7  and  so 
df  df  df 

iffl2  -  k*  -  <S>2  ♦  <3f)a  *  <3f>2  -  #2 

Hence  V®  =  liA  ^  <3l)2  •  (3-U) 

If  A(f)  satisfies  equations  (3.I0)  and  (3.II),  then  (3.O8)  will  be  a  valid 
relationship  on  the  free  streamline.  Since  all  the  functions  involved  are. 
analytic,  this  equation,  valid  for  real  values  of  f,  will  remain  valid  for 
complex  values.  The  flow  given  by  the  resulting  complex  potential  will  have 
the  property  that  the  pressure  is  constant  on  the  streamline  given  by  set¬ 
ting  ip  =  0.  It  will  thus  represent  some  gravity  flow  with  this  as  free 
streamline.  This,  in  essence, is  the  type  of  argument  which  underlies  all 
the  methods  of  this  chapter.  It  depends  heavily  on  the  fact  that  the  free 
streamline  is  necessarily  an  analytic  curve.  This,  of  course,  was  proved 
in  Chapter  Two. 


When  the  free  surface  is  prescribed,  this  method  gives  the  corres¬ 
ponding  flow  under  gravity.  Since  the  free  surface  is  analytic,  there  exists 
a  relation  between  its  slope  and  its  elevation.  Suppose  this  is 
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dx 

dy 


=  cr(-y)  =  or(A)  . 


(3-12) 


Here  cr(A)  is  real  for  real  arguments  and  analytic.  Now,  on  if/  =0,  A  and 
f  are  real  and  7^^  >  (^j)2  by  equation  (3.II).  Hence  equation  (3.O8) 
gives 

x  =  /  E2iA  -  df  • 


Since  on  ip  =  0,  A  =  -y,  this  yields 

I  -  '  'liA 


dA%Pli  df 


Vdf'  J 


dA 


and  so 


’  ^  S  -  Eg iA  -  ^)2)4 


df' 


(3.13) 


Solution  of  this  and  integration  yields 

f.  f[2*(U«*(A)]*dX. 

J 

Substitution  of  equation  (3.I3)  into  equation  (3.O8)  gives 


(3.14) 


z  =  -  iA  -  /  cr(A)dA 


(3.15) 


Equations  (3*14)  and  (3.15)  give  f  and  z  in  terms  of  the  parameter  A 
when  a( A)  is  known.  Again  they  are  valid  for  complex  A  since  all  the 
functions  concerned  are  analytic  and  they  are  consistent  with  Bernoulli 9 s 
equation  on  if/  =  0.  They  provide  the  most  direct  method  known  to  the 
author  of  obtaining  the  flow  having  a  particular  free  streamline,  and  will 
be  re-derived  in  a  slightly  more  general  form  in  the  next  chapter.  The 
extent  of  the  motion  so  defined  depends  on  the  properties  of  cr(A)  for 
complex  values  of  A.  To  quote  Lewy:" [the  motion]  cannot  be  expected  to 
depend  continuously  on  the  shape  of  the  free  surface  since  an  approximation 
to  cr(A)  for  real  A  by  other  analytic  functions  need  not  be  an  approxim¬ 
ation  for  complex  values  of  A,  and  it  is  these  that  determine  the  motion 
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below  (or  above)  the  surface."  This  remark  will  be  validated  somewhat  in 
Chapter  Five  by  some  examples  all  of  which  display  a  singularity  in  the  flow. 


Equation  (3-14)  can  be  re-written  to  give  f  directly  as  a  function 
of  z  provided  A(z)  can  be  determined  by  i.uv  equation  (3. 15). 

Differentiation  of  equation  ( 3 . 1 5 )  gives 


Hence 


[l+cr2]  dA 


=  [f ^T2]^  dz 
L ( -l-cr)2 

=  trfrf^ 

=  U+2i  £]*  dz 


This  gives,  using  equation  (3. 14),  that 

f  -  J  [2gA(  1+2 i  r~)]2  dz  .  (3.16) 

This  gives  f  in  terms  of  z  if  A(z)  is  known. 


Section  3.2  The  Method  of  John. 

In  a  paper  in  1953>  Fritz  John  [3]  outlined  a  method  of  obtaining 
non-steady  gravity  flows.  We  will  restrict  ourselves  to  the  case  of  steady 
flows.  His  approach,  like  that  of  Lewy,  consists  of  writing  down  analytic 
parametric  expressions  for  f  and  z  which  are  valid  on  the  free  stream¬ 
line.  He  chooses  as  parameter  the  time  taken  by  a  particle  to  travel  down 
the  free  streamline  from  some  fixed  point  to  its  present  location.  This 
parameter  is  real  on  the  free  streamline  and  John  analytically  extends  f 
and  z  into  the  fluid  by  letting  it  assume  complex  values. 


Suppose  then  we  are  given  real  parametric  equations  for  the  free 


"  ;V  :p.»;  /  ■  ,  ;;  . 


‘ 


streamline,  the  parameter  being  the  time  t.  That  is,  suppose  xft)  and 
y(t)  give  the  position  of  a  particle  on  the  free  streamline  at  time  t. 
These  are  analytic  functions  of  t  as  was  shown  at  the  end  of  Chapter  Two, 
and  we  suppose  that  they  are  real  for  real  arguments.  If  we  define 

z(t)  =  x(t)  +  iy(t)  (3- IT) 


then  z(t)  gives  the  free  streamline  as  an  analytic  curve  in  the  complex 
plane.  We  assume  t  is  a  real  quantity  on  ip  =  0.  Now: 

i  =  -pi* 

=  ™  exp[ i  tan"1  (~^) ]  , 

where  ds  =  \/ dx^+dy2  is  an  increment  of  arc  length  measured  along  the  free 

ds 

streamline.  Hence  ”  is  the  speed  of  the  particle  along  the  free  stream- 

QL 

line  and,  if  we  denote  the  x-  and  y-  components  of  the  velocity  of  the 
particle  by  u  and  v  respectively,  we  have  on  the  free  streamline: 

=  u  +  iv  .  ip  =  0  ( 18) 

We  will  not  call  this  the  complex  velocity  of  the  particle  but  will  reserve 

df 

this  title  for  the  quantity  oo  =  —  =  u-iv.  It  can  be  seen  in  a  similar 
way  that  the  acceleration  of  a  particle  at  the  point  z(t)  on  the  free 
streamline  is  given  by  ~~ 


Using  this  formulation  of  the  problem  we  can  now  obtain  expressions 
for  z(t)  and  the  complex  potential  f(t)  inside  the  fluid.  We  first 
derive  a  differential  equation  for  z(t.)  by  considering  the  equation  of 
motion  of  a  particle  at  the  point  z(t.)  on  the  free  streamline.  It  is 
acted  upon  by  two  external  forces:  gravity  and  the  pressure  of  neighboring 
particles.  If  gravity  acts  in  the  negative  y-dlrection  the  equation  of 


r  ■ 


/  !» 


,  '  ■*,;.*  •  *->  ;  "...  ■  ■;  . 
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motion  can  be  written  as 


dfz 

dt2 


=  -ig  + 


where  SJ  P  represents  the  pressure  gradient  and  p  is  the  density  of  the 

fluid.  Now  since  the  pressure  is  constant  along  the  free  streamline,  VP 

dz 

is  normal  to  it.  We  also  have  that  i  —  is  a  normal  vector  to  the  free 

at 

dz 

streamline  since  —  represents  the  velocity  there.  It  follows  that  VP 
dz 

and  i  —  are  proportional;  the  proportionality  factor  must  be  real  but 


need  not  be  constant  and  could  vary  with  time.  Call  it 


dt 


(X(O).  Then; 


d2z  .  ,  dv  dz 

dt5  +  1S  =  1  dt  dt 


i/  =  o 


The  function  must  be  real  for  real  arguments  so  as  to  preserve  the 

vector  character  of  equation  (3.I9)  on  the  free  streamline  (where  t  is 
real).  The  equation  can  be  solved  to  yield 


:(t)  =  -ig  f 


ix(*) 


e"1*^  dp 


dh, 


{ 3  <■  20  ) 


If  x(fc)  is  given,  equation  (3.20)  gives  the  free  streamline  for  some 
gravity  flow. 

The  complex  potential  f(t)  can  now  be  determined  from  z(t). 
The  complex  velocity  is 


If  ■  u- iv  ■  z(t)) 

on  the  free  streamline,  using  equation  (3. 18).  Since  t  is  real  on  if/  =  0, 
t  =  t  and  this  can  be  written 


1  =  3^)]  • 


Tp  =  0 


It  follows  from  the  analyticity  of  z(t)  that  —  [z(t)J  will  be  an  anal¬ 
ytic  function  of  t  even  if  t  is  allowed  to  assume  complex  values. 


»V 

li  '  :  7  .  g 
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_  n  .  df  df  dt  df  /dzx-1 

Finally,  since  —  =  —  —  =  —  (— )  ,  we  have 

dz  dt  dz  dt  'dt' 

ft  -fcWtJlSt  W*)]  •  <5-21> 

This  is  a  valid  statement  with  t  real  and  since  z(t)  is  analytic  it 
remains  valid  when  t  assumes  complex  values.  Hence  if  we  define  f(t)  by 
(3-21)  for  complex  t  the  result  will  be  analytic  in  t  and  hence  in  z 
and  will  thus  represent  the  complex  potential  for  some  gravity  flow  consistent 
with  the  free  streamline  given  by  z(t). 

John  uses  equation  (3.I9)  to  show  that  if 

y  =  n(x)  (3.22) 


is  the  equation  of  the  free  streamline,  then  x  can  be  obtained  as  a  function 
of  t  by  solving  a  second  order  ordinary  differential  equation.  On  ip  =  0, 
x(t)  and  y(t)  are  real  and  so  equating  real  and  imaginary  parts  of  equa¬ 
tion  (3.I9)  leads  to  the  equations 


d£x  _  d^  dx 
dts  ~  dt  dt 


and 


&X  dx 
dt  dt 


If  we  eliminate 


dt 


and  write  0  -  g  &♦<£)«  0 


M  ,  /illNSidfx  dri  dfn  ,dxv2  ±1 
11  +  W  Jdt2  +  dx  dP  W  +  g  dx 


then  we  obtain 

o.  (3.23) 


This  is  non-linear  except  where  the  free  streamline  is  a  straight  line. 

Section  3.3  Tangent  Flows  of  Milne -Thomson. 

In  1959  L.  M.  Milne-Thomson  [b]  gave  the  following  extension  of  the 
method  of  John.  He  obtains  a  gravity  flow  which  reduces  to  a  known  gravity- 
free  flow  if  g  approaches  zero. 


K 

-  -  '  '•■'■•OtJ  :■  ;  j  il  :j,&7  1  i ',:w 

i 1  m  ",  ■■  -r", 


-  •'  ■r1  '  '  ,  V  T '  '  ■  .  ;■*  v:  :  ;  .c-;-.-.  y  .  ■. 
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Suppose  we  are  given  the  real  functions  x^(t)  anc*  y^(t)  of 
time  which  give  the  position  of  a  prticle  on  the  free  streamline  of  a  gravity- 
free  flow  at  time  t.  If  we  write 


«0(O  =  xQ(t)  +  i  yQ(t)  , 


(3.24) 


then  Bernoulli's  equation  on  ijj  =  0  reads 


dx^  dy^  dz.  dz, 

^  dt'  +  Mt  '  "  dt  dt 


‘  dt 


=  u: 


(3-25) 


since  g  =  0.  is  a  constant  which  represents  the  velocity  of  the  fluid 

on  the  free  streamline  of  the  gravity-free  flow.  As  we  have  defined  it, 
z^(t)  must  satisfy  equation  (3-19)  with  g  =  0.  That  is 


<1%  _  i  ^0 

-  “  dt  dt 


3.26) 


dt 


dz. 


This  gives  ^  =  -i  ~  [in(^^)].  Since  equation  (3.25)  implies  that  V 

is  the  fluid  velocity  on  if/  =  0  in  the  zero-gravity  case,  we  can  write 
dz  i0. 

d r  =  uoe  and  80 


^  =  -i  — 
dt  dt 


d0  dy 

,  +10  ]  =  ~r~~  =  “T  tav  '  ('~T~) 

0  0  dt  dt  Mx.  ' 


This  is  real  when  t  is  real  and  so  it  follows  that  this  value  of 
be  used  in  equation  (3-19)  to  give  a  gravity  flow.  In  fact 

d2z. 


d2z 

dP 


+  l 


0 


dt- 


dz^  , 
_ O-i  dz 

dt  J  dt 


<Lx 

dt 


can 


(3.2?) 


gives  a  gravity  flow  which  reduces  to  z(t)  =  Az^(t)  +  B  when  g  =  0.  More¬ 
over  this  flow  contains  an  extra  arbitrary  constant  -  the  skin  speed  B  of 


the  gravity-free  case.  This  flow  is  called  the  tangent  flow  to  z  (t). 


If  we  use  equation  (3.25) >  then  equation  (3.27)  can  be  written 


' '  '  '  -  ■  1  ■' i  :  k  '  i  r»j 

, 
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d  rdz 
dt  Ldt 


dz^  dz 

_ 0-.  _ 0 

dt  '  2  dt 

0 


(3.28) 


This  has  the  property  that  it  has  the  same  form  if  another  parameter  t  is 

used  in  place  of  t.  If  zq(t)  =  xq(t)  +  i  Yq(t)  anc^  ^  T  rea^  on  t^ie 

free  streamline,  then  equation  (3.28)  becomes 

j  j  /  dz^  .  dz 

±-  fdz  /  __0  ,  Z2&.  __Q 

dr  Ldr  /  dr  ^2  dr 


This  can  be  solved  in  the  form 


z  =  A  z  +  B  -  r 
0  u 


0 


,  z  dz^  , 

J  0  0 


(3.29) 


where  A  and  B  are  constants  of  integration.  The  speed  V  on  the  free 
streamline  of  the  gravity  flow  is  obtained  as  follows'  From  equation  (3< 2 
we  have 

3f-£  • 

Using  equation  (3.25)  it  follows  that 

=  “US  l&!S  -  Xm{AZo)  +  §1  U0I2.  (3-30) 

Equation  (3.29)  yields  the  following  expressions  for  x  and  y.  If 


B  =  b^  +  ib^,  we  have 


X  -  Re(Az0)  +  +  5§-  /  (x0dy0  - 


and 


y  =  Im  (Az  )  +  b 


— I-  Iz  I2 
O'  r  “2  2U2q  1  01 


The  second  of  these  and  equation  (3*30)  imply 


V2  +  2gy  =  Ug  | A | 2  +  2gb2  . 


KT&'S  Ul  O-tV:  &<$  fifiO 
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This  has  the  form  of  equation  (2.01)  if  we  choose  the  origin  of  coordinates 

uoH2 

such  that  b  =  -  — - -  .  We  can  choose  b  =  0  and  so  get 

1 


and 


x=Re(Az0)  +  ^  /<x 0  %  '  y0  V  (5‘3l) 

Un  lA|2 

y  =  Im<Azo>  -  iu®  l20lS  -  — zT  ■  (3-32) 


If  we  choose  the  free  streamline  in  the  gravity-free  case  to  be 
a  circle,  then  equations  (3*31)  and  (3-32)  give  the  example  done  in  John’s 
paper  [3]  where  the  free  streamline  is  a  trochoid.  In  addition  spiral 
shaped  curves  can  lead  to  monotonic  tangent  free  streamlines.  In  fact  the 
shape  of  the  tangent  free  streamline  seems  in  general  to  bear  little  resemb¬ 
lance  to  that  of  the  free  streamline  in  the  gravity  free  case.  For  this 
reason  it  would  seem  that  the  name  Tangent  Flow  is  a  rather  inept  one. 
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CHAPTER  FOUR 

RELATIONS  BETWEEN  THESE  FORMULATIONS 

It  is  evident  that  the  above  are  all  different  formulations  of  the 
solution  to  the  same  problem.  To  indicate  precisely  the  connection  between 
them  we  first  reformulate  the  problem  and  obtain  Lewy's  results  in  a  slightly 
more  general  form,  and  then  derive  the  connections  between  these  and  the 
results  of  John  and  Milne -Thomson. 

If  we  choose  the  quantity  A  such  that 

A  =  -y  7p  =  0  (4.01) 

on  the  free  streamline,  then  the  Bernoulli  equation  can  be  written 

V2  =  U2  +  2gA  (4.02) 

there,  where  U  is  a  constant  with  the  dimension  of  velocity.  Lewy  found 

it  convenient  to  choose  it  equal  to  zero  but  we  leave  it  arbitrary.  Equation 

U2 

(4.02)  requires  that  A  >  -  — —  on  the  free  streamline. 

We  assume  that  we  are  given  a  real  analytic  curve  which  we  want 
to  be  the  free  streamline  of  a  gravity  flow.  Suppose  it  has  the  form 

x  =  x(A)  =  x(-y)  .  (4.03) 

If  we  use  equation  (4.0l)  this  curve  is  given  by 

z  =  -iA  +  x(A)  (4.04) 

in  the  plane  of  the  complex  variable  z  =  x  +  iy.  As  in  Lewy’s  paper  x(A) 

is  real  for  real  arguments,  and  so  equation  (4.04)  represents  the  free  stream 

U2 

line  if  A  is  real  and  A  >  -  —  . 


■  . 
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We  want  to  construct  an  analytic  velocity  function  u>  with  the 
property  that  w  to  =  V2  =  U2  +  2gA  on  the  free  streamline  in  accordance  with 
equation  (4.02).  Since  A  is  to  be  real  on  the  free  streamline,  the  function 


w  =  (U2  +  2gA)2 


satisfies  this  condition  if  0(A)  is  an  analytic  function  which  is  real  for 


U2 

real  arguments  A  >  -  —  . 


Now  suppose  f  is  the  complex  potential.  Then  since  w  = 
we  have  using  equations  (4.04)  and  (4.05)  that 


df  df  dz  ,  .  dxv 

-  “(-1  +  • 


dA  dz  dA 


df 

dz 


Hence: 


f  = 


(U2+2gA)2  (cos0+i  sin0)  (-i+^)dA  =  cp(A).f 


Now  if  we  call  the  free  streamline  ^(A)  =  0,  then  since  A,  @(A)  and 
x(A)  are  all  real  there,  we  can  write 


^(A)  =  f  (U2+2gA)2  (sin0  ^  -  cos@)dA  =  0 


J 


on  the  free  streamline.  This  requires 


dx  _  -l,dx%  . 

Z  =  cot(®)  »  or  Q  =  cot  fe) 


.06) 


dA  -  — -v-/  >  ^  -  -  —  vdA' t 

This  clearly  gives  w  the  direction  of  the  tangent  to  our  curve  which  is 
as  it  should  be. 


Following  Lewy  we  now  define  the  analytic  function 


*<A )--£  • 


(4.07) 


Again  cr(A)  is  real  if  A  is  real  and  A  > 
we  can  write  equation  (4.04)  as 


uf 

2g 


.  Using  equation  (4.07) 


■  v  '  "i;>:  ■:  X  »  *  . :  f  yx  . ;  i  sy-  ■„  yo 

v  I  r 
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z 


=  -iA  - 


r  ci(A)dA 


which  is  equation  (3*15)  of  Lewy.  We  have  also 


i@(A)  i  cot  (-O')  ,0-1 

e  =  e  '  =  ( — -7 

vcr+i 


Hence: 

“  -  (U2+2gA)=  (2=i) 2  . 

P  1  1 

and 

f  =  /  (U2+2gA)2  (l+cr2)2dA 

(4.08) 


(4.09) 

(4.10) 

(4.11) 


Equation  (4.11)  reduces  to  equation  (3.14)  of  Lewy  if  U  =  0  and  so  pro¬ 
vides  a  slight  generalization  of  that  equation.  The  factor  U2  in  equations 
(4.10)  and  (4.11)  allows  us  to  construct  flows  which  have  the  property  that 

they  are  uniform  streams  of  finite  depth  at  great  distances  upstream.  In 

1 

addition  the  factor  (cr+i)2  in  the  denominator  of  equation  (4.10)  shows 
that  a  singularity  in  the  flow  is  likely  to  occur  where  cr  =  -i„  This  branch 
point  in  the  velocity  determines  a  natural  bound  on  the  extent  to  which  the 
flow  can  be  analytically  continued  away  from  the  free  streamline,  a  phenomenon 
which  is  in  keeping  with  the  remark  of  Lewy's  quoted  in  Chapter  Three,  Some 
examples  of  this  will  be  3een  in  Chapter  Five. 


We  next  obtain  expressions  for  the  quantities  t  and  x  °f  John 
in  terms  of  A  and  cr(A).  We  begin  by  obtaining  equation  (3. 21)  using  a 
method  different  from  that  used  by  John. 


form 


If  we  write  w  =  u  -  iv,  then  we  can  write  equation  (4.02)  in  the 


U2  +  2gA  =  w  w  =  u)(u  +  iv) 

Following  John  we  define  a  new  parameter,  the  time  t,  on  the  free  streamline 
by 
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dz 

dt 


=  u  +  iv  =  w  . 


(4.12) 


We  thus  have  co(u+iv)  =  4^  •  ~ 

'  dz  dt 

free  streamline  becomes 


~  and  so  the  Bernoulli  relation  on  the 
dt 


f  +  2gA 


tf/  =  0 


13) 


This  implies  that  the  time  t  is  a  real  quantity  on  the  free  streamline  since 
f  and  A  are  real  there.  Using  these  results  we  can  write 


df  -  dz  dz 

dt  -  w  w  =  dt  dt 


on  the  free  streamline.  Since  t  is  real  there,  then  if  we  refer  z  to 

dz  _  d_ 
dt  "  dt 

(4.14) 


an  argument  t  we  can  write  z(t)  =  z(t)  and  so  =  ~r  z(t).  Finally 


df  d  d  fZ\ 

z(t)  ^  *(t)  , 


dt  dt 


which  is  equation  (3.2I)  of  John. 

The  connection  between  the  quantities  A  and  t  on  the  free 
streamline  is  found  from  equations  (4.11)  and  (4.I3).  We  have 

3F  -  U2  +  2g?v  =  ff  .  (Us  +  2g*)T  (1  +  O*  57  • 


Hence : 


C  d?V 

which  is  an  analytic  function  of  A  on  the  free  streamline. 


-.15) 


We  are  now  in  a  position  to  write  down  x  as  a  function  of  A 
and  O'(A).  John  writes  the  equation  of  motion  of  a  particle  on  the  free 
streamline  as 

+  ig  =  i  ^  — 
dt2  +  g  1  dt  dt 


' ■;  ■ 

. 


i  r  I  SCO  :  ;!|  '  v  .  - 
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where  the  function  x  is  real  on  the  free  streamline.  We  can  write 

d2z  d  ,dz,  dz  -  dw 
dt2  ~  dz  dt'  dt  “  dz 

using  equation  (4.12),  and  so  the  equation  of  motion  becomes 

-  du>  .  ,  dv  ,-.2 

U  dl  +  1S  =  1  dt  ('0)  * 

Some  algebra  and  integration  gives 

X  =  ~i  •Mw)  +  T ||y2  . 

Now  from  equation  (4.08),  dz  =  (“i~cr)dA  and  so,  using  equation  (4.05), 

X  -  -i[-i0+Jin(U2+2gA)]  +  P  dy  . 

Using  equation  (4.09)  we  obtain: 


=  -e-  “  dn(U2+2gA)  -  P  (~r)dA 

^  J  U  -i-clgA  oq*i 


Finally, 


v  „  „  p>  ^  P 

X  "  J  U2f-2gA 


Xl6) 


Using  equations  (4.06)  and  (4.07)  we  -an  write  this  as 


-i,  %  r  g^A 

X  =  cot  (cr.  -  /  T^-TT’T 
'■  J  U*--f2gA 


(4. 1?) 


This  again  is  an  analytic  function  of  \  in  the  free  streamline, 


We  now  establish  an 

Thomson  and  the  quantities  A 

of  time,  z^(t),  satisfying 

dz^  dz^ 
_ 0 _ 0 

dt  dt 


analytic  relationship  between  z^  of  Milne- 
and  x*  If  we  are  given  an  analytic  function 


=  U2 


(4.18) 


then 


’  ■ r.)  .  n.s~  £/• 


-  27  - 


=  Az  +  B  - 

0  U2 


z„  dz 
0  0 


(4.19) 


is  a  gravity  flow.  Equation  (4.18)  is  just  equation  (3.25)  with  chosen 
equal  to  U  as  defined  by  equation  (4.0l).  Differentiating  and  using  equa¬ 
tion  (4.12): 


dz 

dt 


"  (A  ■  zo}  3T  "  5 


Hence: 


U) 


- (i  +  u§  Vir  = uS^ +  u§  z0y 


dt 


using  equation  (4.18).  This  gives 

dz. 


0 


dt 


w  U2A+igz 


0 


If  we  use  equation  (4.13),  then 


df  df  dt  ,tt2  _  -wdt 

“  =  33  =  33  57  =  (u  +  2sX)3I 


and  so 


dt 

10 


dz 


dz. 


0 


U2+2gA  U2A+igz 


Since  dz  =  (-i-cr)dA  from  equation  (4.08),  then  we  can  write  this  as 


dz. 


_____  -  ( i-t-cr)  dA 

1  +  z0  1  +  P  * 


where  the  denominators  are  dimensionless.  Integration  gives: 


ii  +  u§  z0)  =  ■  fir  <n(1  +  P  ■ 


cr  d  A 

1  +  p  A 


Jtn 


ft  +  pg  z0  1  r  sad* 

1  ♦  §§  *)41  =  j 


This  reduces  to 


- 
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If  we  use  equation  (4.16)  to  eliminate  the  integral  on  the  right  we  obtain: 


i  +  =  <1+§#> el(x+e) 


Finally,  using  equation  (4.05),  we  have 


1  +  u§  z0  =  n  elx 


,20, 


We  can  now  show  that  the  constant  A  in  equation  (4.19)  can  be 
chosen  to  be  zero  without  loss  of  generality.  This  merely  amounts  to  choos¬ 
ing  a  different  function  z^.  To  see  this  suppose  we  choose  a  function  z^ 


given  by 


i§  z  -  A  +  iS  z 
U3  1  ~  +  U2  0  * 


dz, 

We  have  that  — -  =  and  so  equation  (4.18)  is  satisfied  by  .  Also 


we  have 


/ 


_  p _  iu2 

z„  dz,.  =  /  z,  dz,  -  -  A  z,  . 

0  0  J  1  1  g  1 


Hence,  from  equation  (4.19)  we  have 


+  B  -  M  f  z  dz  -  Az 


1 


-*-si/wT‘ 


This  is  just  equation  (4.19)  with  A  =  0,  z^  =  z^,  and  the  origin  of  z 
translated  vertically.  Hence,  the  simplest  formulation  of  the  tangent  flows 
is: 


Ifo  l!o 

dt  dt 


IT 


2  -  c  -  u§  fzo  dzo' 


iUu)  iv 

z  =  -  -  e  A" 

0  g 


.21) 


where  C  is  a  complex  constant. 
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Equation  (4.20)  shows  that  z^  is  an  analytic  function  of  A 

since  u>  and  x  are.  Also,  every  flow  determines  a  cr(A)  which  is  anal- 
U2 

ytic  for  A  >  -  —  .  This  in  turn  yields  a  z^.  Conversely  every  z^ 
determines  some  flow.  Hence,  every  flow  is  a  tangent  flow. 

Equation  (4.21)  admits  another  conclusion.  If  we  replace  z^ 
by  zq^  where  a  is  a  real  constant,  then  the  tangent  flow  with  free 
streamline  given  by  z  is  unchanged.  Hence  the  orientation  of  a  given 
gravity  free  flow  leaves  the  resulting  tangent  flow  unchanged,  and  so  we 
can  say  that  every  flow  is  tangent  to  an  infinity  of  flows. 

Since  the  function  x  appears  in  a  differentiated  form  in  equa¬ 
tion  (3. 19),  then  two  functions  x^  and  differ  by  a  constant 

will  lead  to  the  same  gravity  flow.  The  function  <j,  however,  has  no  such 
arbitrary  character. 

We  now  cite  an  example  which  has  a  free  surface  consisting  of 

2it 

two  straight  lines  meeting  in  an  angle  of  ~  .  We  will  use  the  solution 
embodied  in  equations  (4.02),  (4.08)  and  (4.11).  However  we  choose  U  =  0 
so  as  to  simplify  the  notation.  This  means  that  the  corner  in  the  free 
streamline  will  occur  at  A  =  0  since  it  is  a  stagnation  point.  The 
rest  of  the  free  surface  will  correspond  to  A  >  0,  that  is  to  ssy  y  <  0. 

If  p  is  a  constant  acute  angle,  and  if  we  choose 


cr  =  -  cot  p  , 

then  equation  (4.08)  with  U  =  0  gives 


z 


( -i  +  cot  p)A  = 


Ae 


-ip 


sinp 


> 


and  equation  (4.1l)  gives 
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£  =  f  (2gX)^  (l+cot2p)^dX  -  A?/2  . 

Eliminating  A  these  give 

313 

f  =  |  (29  sinpp  e  2  z5/2  # 

If  z  =  re*^  ,  then 

f  =  f  r?/2  (2g  sin3)i  e^W) 

The  free  streamline  corresponds  to  tJ>  =  0  and  this  occurs  when 


f  (@  +  3) 


=  n  jt 


9  =  -3  + 


2njc 


Since  y  =  r  sin  9  ,  the  condition  that  y  <  0  implies  that  -n  <  6  <  0, 
Hence  only  two  acceptable  values  of  9  are  possible;  choose  them  to 
correspond  to  n  =  0  and  n  =  -1.  These  give 


M  , ,  df 
Now  c*)  =  — 
dz 


9  =  -3  , 

5i3 

JL  i 

(2g  sinp)2  e  z2 


0  =  -(3  +  “) 


Hence 


V2  =  w  w  =  2g  sin3  | z |  =  2g  sin  3  r 


With  U  =  0,  Bernoulli's  equation  (i+,02)  reads  V2  =  -2gy  on  the  free 
streamline.  In  the  present  example  this  gives 

2g  sin3  r  =  -  2gy  =  -2g  sin  9  r  . 


The  portion  of  the  free  streamline  given  by  9  =  -3  satisfies  this. 
The  portion  given  by  0  =  -(3  +  ~)  does  so  if 

sinS  =  sin(3  + 

3 


■ 

/  ■  ,  -•  j :  y.  ■  ■  ; 
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This  implies 


Hence 


f 


2  nTr 

3 


H  z3/2  =  2j/g  r3/2 
3 


♦  f; 


and  so  the  stream  function  is  given  by 

,  2  \Tg  3/2  .  ,39  it. 

i/  =  -y®  r  '  sin(  ^  +  p  ♦ 

If  we  set  0  =  -  ^  +  ^  >  then  sin(-^  +  ^)  =  -cos  ^  .  Thus 

*  -  -  ^  cos(f)  . 

Tt 

The  two  free  streamlines  correspond  to  £,  =  +  ^  .  The  streamlines  correspond 
to  ip  =  -  ipQ,  where  ip^  >  0.  Their  equation  is,  then, 

r  =  (2/f^2/3  sec2//5  • 

The  highest  point  on  any  of  these  curves  occurs  when  |  =  0;  that  is  when 
9  =  -  ~  .  They  are  symmetrical  about  this  point.  A  sketch  of  the  flow 
appears  in  figure  4.01  with  a  typical  streamline  included. 


Figure  4.01 
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CHAPTER  FIVE 

EXAMPLES 


The  aim  here 
like  a  uniform  flow  of 


is  to  construct  some  examples  of  flows  which  behave 
finite  depth  at  great  distances  upstream. 


We  have  that 


U) 


the  complex  velocity  is  given  by  equation  (4.10)s 

-  (u2  +  t$j  ;■  ^  •  (5.01) 


Since  cr(A)  is  real  for  real  A  >~  ~  ,  then  Jw|  -> U  as  A  0 

through  real  values.  Moreover ,  since  c(A)  =  on  the  free  streamline  we 
must  have  cr(A)  — ►  3-ar  upstream  if  the  streamline  is  to  be  horizontal 

there.  Hence,  if  we  arrange  it  that  A  =  0  corresponds  to  a  point  far 
upstream  on  the  free  streamline,  then  we  will  have  a  uniform  flow  there. 

In  addition,  if  it  can  be  shown  that  the  total  flux  of  fluid  between  some 
particular  streamline  and  the  free  streamline  is  finite,  then  it  follows  by 
continuity  that  the  depth  of  flow  far  upstream  is  finite. 


If  we  restrict  our  attention  to  those  flows  for  which  the  free 
streamline  is  continually  falling,  then  cr(A)  <  0  for  all  A  >  0  .  In 

particular  we  consider  flows  which  have  o*(A)  — >  ~  oo  as  A  -+  0,  The 

streamline  is  given  in  the  complex  z-plane  by  equation  (4„08)s 


■iA 


dA 


(5.02) 


In  most  of  our  examples  we  will  consider  flows  for  which  Jcr{l\)  dA  -►00 
as  A  -*  0.  The  form  of  the  streamline  in  the  z-plane  is  shown  In  the 


sketch  of  figure  5.01. 


'  -  ■’  '  '  .  >'■  •  ..  ■■  '  •:  .  ■■■  ;  V  ,  •  j. 

■ 


■  ,  '  4*6  -1 

■  VrC.  ■  .  V  '  .  :  :  i'i  v  ,.-M 

!"  ■  '  "  '■ . "  -  .  :  V  ■  :  '  O  tB 
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Figure  5.01 


It  will  now  be  shown  that  values  of  z(A)  lying  beneath  the 
free  streamline  correspond  to  Im(A)  <  0.  We  first  fix  the  position 
of  the  y-axis  along  the  free  streamline  by  writing  equation  (5«02)  as 


:(A)  =  -iA  -  J  ct(A)  dA  , 


:5.03) 


where  c  is  real  and  positive.  Now,  if  we  are  given  an  arbitrary,  real, 
positive  number  a  ,  then  since  the  free  streamline  is  analytic,  cr( 
is  analytic  near  A  =  a  .  Moreover,  cr(a)  <  0  .  Hence, 


<r(A)  =  cr(a)  +  o-'(a)  (A-a)  +  icr"(a)  (A-a)2  +  ... 


for  A  sufficiently  near  a  .  Equation  (5.03)  then  gives 


z(A)  =  -iA  -  f  ct(A)  dA 
J  c 

r  a  r  A 

=  -i(A-a)  -  ia  -  /  ct( A)  dA  -  /  o-(A) 

J  c  J  c 

r a  ip 

=  -ia  -  J  cr(A)  dA  -  i(A-a)  -  [cr(a)(A-a)  +  2cr'(a)(A-a)  +  ...] 


dA 


z(a)  -  (i  +cr(a))(A-a)  -  20-'  (a)(A-a)2  +  ... 
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Now  suppose  A  =  a  +  ie  where  €  is  small  and  real.  This  gives 

i  2 

z(a+ie)  =  z(a)  +  (l-icr(a))e+  2ct' (a)e  +  ... 

Since  a(a)  <  0  ,  this  shows  that  for  small  negative  €  ,  z(a+ie) 
lies  below  and  to  the  left  of  the  point  z(a)  on  the  free  streamline. 
That  is  to  say,  z(a+ie)  is  in  the  fluid  below  the  free  streamline  if 
€  =  Im(A)  <  0  , 


From  equation  (5. 01)  it  is  evident  that  if  <t('a)  =  -i  ,  then 
there  will  be  a  singularity  in  the  velocity  provided  U2  +  2gA  j  0 
at  the  same  point.  Define  and  z by* 


o-(^c)  =  _i  i 


and 


=  z(A  ) 
c ' 


|  (  5.04  ) 


If  we  assume  that  cr(7\ )  is  analytic  near  Ac  ,  then  we  can  determine  the 
nature  of  the  singularity  there.  Since  0(A)  =  -i  ,  we  can  write 


a(\)  =  -i  +  cr'(7\c)(A-Ac)  +  .. 


for  values  of  A  sufficiently  near  A  ,  We  have  from  equation  (5.03) 

2  C 

d  z  d  z 

that  =  -i  -  cr(A)  and  — =  ~<x'(A)  .  The  first  of  these  gives 

dA 


=  0  when  A  =  A^  .  Equation  (5.02)  implies  that  z(A)  is 

analytic  near  A  =  A  if  O'(A)  is.  Hence 


z(A)  =  z  -  2cr'(A  )(A-A  )2  +  ... 


near  A  =  A  .  We  now  write 
c 


A  =  A  +  pe 
c 


iS 


and  consider  the  situation  when  p  <v  0  .  The  expansions  of  z(A)  and 


' 
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cr(A)  become; 


=  -i  +  cr' (A  )pei&  + 


z  -  z 


i  ,  /  ,  x  2  2i5 
=  -2cr'  (  A )p  e  + 


when  p  0  .  In  this  notation  equation  (5.0l)  becomes 

'-2i  +cr'  (A  Jpe^Nli 


to 


-[<* 


+  2gA  +  2gpei5) 


cr'(A  )pe 


i8 


— -  (U2  +  2gA  +  2gpel5)(l  4-  ^cr!  (A.  )pei6) 


“pei8cr'  (Ac) 


.5  it 

12  2  + 

-  *  p  e 


[-^(h(u£  +  2*c>. 


where  we  are  keeping  only  the  leading  term  on  the  assumption  that  p  is 
small.  The  case  where  U2  +  2gAQ  =  0  will  be  dealt  with  below. 


and 


then 


If  we  now  define 

P  =  arg[-o!(Ac)] 
v  =  argtU2  +  2gAc ]  , 


1 


I  (5.05) 

i(v  -  p)  .  Hence 

-  p)  . 


Also  arg(z-zc)  =  25  +  P 

Since  arg(co)  is  the  direction  of  the  fluid  velocity,  then  the  streamline 
through  zc  will  correspond  to  values  of  6  for  which  the  quantities 
arg((o)  and  arg(z-zc)  are  equal  or  differ  by  a  multiple  of  x  „  To  find 

a 
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these,  set 

25  +  P  =  g  -  £  -  2  ( v  -  p  )  +  njr  . 
That  is  5=-^(v+|i)  +  “g~  Jt  . 


The  three  independent  values  of  5  correspond  to  n  =  0,  1,  2  J  all 
others  differ  from  these  by  a  multiple  of  2jt  .  We  have  that 


and 


arg(z-zc)  =  ^(p  -  2v)  +  it  , 

arg(uj)  =  i(p  -  2v)  +  (^)  *  . 


8ir 

As  n  increases  by  two,  arg(z-z  )  increases  by  ~  and  so  these  two 

c  3 

values  do  not  lie  on  the  same  physical  plane.  Hence  the  streamline 
entering  and  leaving  the  corner  must  correspond  to  two  consecutive  values 
of  n  .  Since  we  want  w  near  the  corner,  and  since  we  want  the 

flow  to  be  from  left  to  right  in  figure  5. 01 >  then  the  pair  n=l,  n=2 
are  inadmissible.  We  have,  then,  for  n  =  0,  1, 


n  =  0  arg(z-zc)  =  arg(w)  =  ^(p  -  2v)  -  |  • 

n  =  1  arg(z-z  )  =  ^(p  -  2v)  +  it  ,  arg(w)  =  ~(p  -  2v)  . 

c  3  5 

The  first  corresponds  to  a  streamline  leaving  z c  ;  the  second  to  one 
entering.  If  we  define  the  angle  c  as 

c  =  |  .  I(p  .  2v)  ,  (5.06) 

then  the  situation  in  the  neighborhood  of  z is  as  appears  in  figure 


5.02. 


■ 


fS 


.  :  ■  -  3  ;  >■[■?!  o  ■  j;j:  3  )n  ob  '  .  .  h 

'  (}  '  ‘  ■  :!,  '  i:\:  :  •  v,  ■-  r  :  t-..;  ..  - :  :i 

■  '  -  :  i  4  .  M 
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Figure  5.02 


The  streamline  enters  the  singularity  and  is  rotated  through  an  angle 

of  -  —  before  it  leaves.  Since  we  cannot  abide  a  singularity  like  this 
3 

inside  the  fluid,  it  is  necessary  to  choose  the  streamline  through  z c 
or  one  above  it  as  the  bed  of  our  flow.  If  we  choose  the  one  through  z 


the  situation  is  that  of  a  flow  over  a  sharp  corner  of  angle 


2k 


In 


the  examples  that  follow  we  will  refer  to  this  streamline  as  the  bed  and 
to  the  point  z as  the  corner. 

2jt 

Corners  with  angles  other  than  —  are  also  possible.  It  is 

3 

the  behaviour  of  <t(A)  near  A  =  Ac  that  affects  the  situation.  If, 
in  fact,  we  suppose  that  near  A  =  Ac ,  cr(A)  behaves  like 


cr(A)  =  -i  +  a(A-Ac)7  , 


(5.07) 


then  other  angles  are  obtained.  Here  a  and  y  are  constants  and 
7  >  0  .  We  have  in  general  from  equation  (5. 03)  that 


r  c 

:(A)  =  -iA  -  /  or(A)  dA 

J  c 

-  J  c  cr(A )  dA  -  i(A-Ac)  -J  <r(A)  dA  . 


=  -iA. 


If  we  use  equations  (5.04)  and  (5.O7)  we  can  write  this  as 
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f  A 

:(A)  =  z.  -  i(A-Ac)  -J  [-i  -  a(A~Ac)7] 


dA 


=  z  -  — 
c  y 


TT  <™c>r+1  • 


Now,  as  above,  if  we  write 


i  6 


A  —  Ac  +  pe  , 


then  we  obtain 


cr(A)  =  -i  +  ap7e17&  , 


and 


zM  =  Zc  "  7  +  1 


—  p7+1  e1^74,1^6 


when  p  ~  0  .  These  relations  together  with  equation  (5.0l)  yield 


oj 


-  [(U2  ♦  2gAc  +  Sgp.18) 


7  iy& 
ap'  e  ' 


(5.08) 


We  now  have  two  situations  to  consider  according  as  the  quantity 
+  2gAc  is,  or  is  not,  equal  to  zero.  The  first  possibility  can  only 
occur  if  Ac  is  a  negative  real  number  and,  at  the  same  time,  the 
velocity  U  far  upstream  has  an  appropriate  value.  One  of  the  examples 
below  exhibits  this  property.  We  will  discuss  both  possibilities. 


Fi 


rst  assume  that  +  2gA£  j  0  . 

rCu2  +  2gX  )(-2ih2 


Then  equation  (5.O8)  gives 


u>  = 


7  iy6 
ap  e 


where  only  the  largest  term  has  been  retained.  If  we  now  define 

Hj.  =  arg(-a)  (5.09) 

and  keep  v  =  arg(U^  +  2gA  )  as  in  equation  (5.05)  ,  then 


« 

■ . ~  s 

■a  +  :f 

<  -f  A  *  A 

■■  '/  0  ,  *  i. 


S  -  ,i  ■  *  '•  }"| 
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arg(ui)  =  p-  -  J  -  (v  -  (i1)  . 

In  addition 

arg(z-zc)  =  (7+1)5  +  . 

If  we  again  reason  that  values  of  5  for  which  arg(to)  and  arg(z-z^) 
differ  by  a  multiple  of  jt  correspond  to  streamlines  through  z^  ,  then 
these  values  are  given  by 

(7+1)5  +  =  2  S  "  I  "  ^(v  -  ^)  +  nr£  i 

^gr  5  =  -  2(v  +  nx)  +  (~f~)  it  , 

v  +  Mi  te  -  1 

5  M  '  TTz~  +  sUW)x  ■ 

This  gives 

arg(z-zc)  =  -  ( v  +  M’1  )  +  | (4n-l)  +  ^ 

=  7+2  ^1  ”  +  |  ^  (4n-l )  . 

If  n  increases  by  two,  then  arg(z-zc)  increases  by  an  amount 

(2-™)(4jt)  which  is  greater  than  2it  for  all  7  >  0  „  Hence  n=0  and 
7+<- 

n=l  give  two  values  of  arg(z-zc)  which  lie  in  the  same  physical  plane. 

The  first  means  arg(z-zc)  =  arg(w)  and  so  corresponds  to  a  streamline 

leaving  z ^  5  the  second  corresponds  to  a  streamline  entering  z^  ,  These 

two  values  of  arg(z-z  )  differ  by  (2jt)(~ — — )  with  the  streamline  leaving 

c  7  +2 

the  corner  corresponding  to  the  smaller  value.  Hence  the  angle  of  the 
corner  itself  is 

c0  =  2*  -  &(£—)  =  .  (5.10) 

This  varies  from  n  down  to  zero  as  7  runs  from  0  to  °o  o  The 


situation  is  as  is  shown  in  figure  5.03. 


-  40  - 


Figure  5.03 


The  angle  c  in  this  case  is 


c  =  2-  (2±1) 

2  >+2; 


-  53^1  -  ^Dv] 


(5.H). 


If  cr(A)  is  an  analytic  function  of  A  near  A  =  Aq  ,  then 

7  will  be  an  integer.  This  being  the  case,  then  if  we  suppose  that 

cr(Ac)  =  -i  ;  and  that  the  r^  derivatives  of  ct(A)  is  the  first  one 

which  does  not  vanish  at  A  =  Ac  ,  then  equation  (5.07)  implies  that 

7  =  r  and  a  =  — -  oA  '(A  )  .  The  computation  leading  to  equation 
^  ^ 

(5.06)  is  the  case  7  =  1.  Equation  (5.H)  reduces  to  equation  (5=06) 

2it 


in  this  case  and  equation  (5.  10)  gives  c  = 

0  5 


All  the  examples  in 
this  chapter  involve  this  case;  that  is  to  say  the  choices  of  the 
function  cr(\)  which  we  will  make  will  be  analytic  at  A  =  A£  and  will 
have  cr  ’  (A  )  ^  0  . 


s' 

A  different  situation  occurs  if  IP  +  2gA  =  0  .  In  this  case 


equation  (5.08)  gives 


U)  = 


i&  \  I  -2i  +  aprei,'5N'"12 
(2gpe  )  (  -  - 

ap  e 


-  '  1  ■'  .  : ;  '  •  ■  :  .  .  "  ...  :■ ■  ■-  ....  <  ;  c  ,?) 

J  .■  •  u  : <n  '  aj>  ?  wj 


0  k  ■'  )  '  t>  •  ' 


0.' 
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If  p  ;v  0  this  is 

to  -  2/=S:  e-i(r-i)s  +  *5  [JLji 

Jp?-1 

which,  taken  with  equation  (5.09)>  yields 

p. 

arg(co)  =  (7-1)5  -  J  +  2^  • 

Again  arg(z-zc)  =  (7+1)5  +  . 

As  before  we  set 

^1 

(7+l)5  +  P-^  =  (7-1)5  -  ]£  +  2~  + 

to  obtain  values  of  5  corresponding  to  streamlines  through  z ^  .  This 
gives 

4n-l 

8  =  '  5“  +  ~ T  *  ’ 

and  so 

arg(z-zc)  =  +  2^1  (4n-l  )jt 

gives  the  orientation  of  the  possible  streamlines.  Again  n=0  corresponds 
to  a  streamline  leaving  z c  and  n=l  corresponds  to  one  entering.  The 
angle  c  is  given  by 

c  =  n  ’  ^1  (5.12) 

Since  arg(z-zc)  increases  by  jt  as  n  increases  by  one,  then  the 

angle  c^  of  figure  5.03  Is  given  by 

c0  =  2*  -  ^  ^  .  (5. 13) 

This  gives  that  c^  ^  «  according  as  7  ]>  1  .  If  7  >  1  we  have 


the  case  of  a  corner  like  those  which  we  have  been  discussing.  If 


[  -  t  C:.;i 

BW  S'JC 

'  r‘s  ■  ;  ..  ic  a*  '•  ,  . 

.  '■  ■: :  ' !■.  .  ?  ,0  V\' 


■  •'  r:  < j!X' 
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y  =  1,  Cq  =  it  and  there  is  no  corner.  Finally  if  7  <  1  we  have 

w  =  0  at  the  corner  and  c  >  jt  .  The  extreme  values  are  y  -  0  and 

7=3.  1°  the  first  case  c^  =  ^  and  so  the  fluid  flows  inside  a 

corner  of  angle  ^  .  If  7=3  then  c^  =  0  and  so  this  is  the  case 

of  flows  over  a  sharp  cusp.  Values  of  7  >  3  give  c^  <  0  and  so 
are  non  physical.  It  is  of  some  interest  to  point  out  that  7  =  2 
gives  Cq  =  ^  independently  of  whether  U2  +  2gA  is  zero  or  not. 


Since  the  velocity  far  upstream  is  known,  then  we  can  calculate 
the  depth  there  if  we  know  the  total  flux  of  fluid.  We  calculate  this 
by  determining  the  change  in  the  stream  function  ip  between  the  bed  and 
the  free  surface.  There  are  two  quite  general  approaches  to  this  problem. 
If  it  happens  that  the  free  streamline  assumes  a  slope  of  -1  at  some 
point  then  cr  =  -1  there  and  so,  since  the  corner  corresponds  to  cr  =  -i, 
we  can  find  the  change  in  f  between  these  points  by  integrating  its 
derivative  along  a  portion  of  the  unit  circle  in  the  cr-plane.  Since,  by 
equation  (4.1l), 


f  = 


(U2  +  2gA)2  (1  +  cr2)2  dA  , 


then  if  we  denote  the  change  in  f  between  the  bed  and  the  free  streamline 
by  Af  ,  we  get 

-1 


Af  =  J  (U2  +  2gA)2  (1  +  cr2)2  dcr 


wherever  we  can  invert  cr  =  cr(A)  to  obtain  A  in  terms  of  cr  .  If  we 
write 


cr  =  -e 


IQ 


r.  2  r  219 

vl  +  cr  =  vl  +  e 


1® 

2 


n/2  cos  6  . 


then 


e 


•  ?  V-  '  •  ,v 

'  T;  -  '  ■  '  .  ::  tt.i.  ■.  :s  ,  l, a...  ; 


■ 
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Hence 


Af 


=  -  I  2  (U2  +  2gA)2  (  cos  0) 


d A 
d0 


d0. 


(5.14) 


0 


The  change  in  the  stream  function,  A^  ,  will  be  the  imaginary  part  of 
this.  Equation  (5.14)  shows  in  addition  that  a  sufficient  condition 
that  the  stream  have  a  finite  depth  far  upstream  is  that 
be  bounded  for  0  <  9  <  ^  where,  of  course,  A  =  A(o")  =  A(-e^)  . 


KtfWo2  j|| 


We  can  also  compute  Af  from  an  integral  in  the  A-plane  in 
many  cases  of  interest.  If  the  corner  corresponds  to  A  =  where 


-ia 


=  R  e 
c 


then,  for  a  wide  class  of  flows,  A  =  Rc  will  correspond  to  some  point 

on  the  free  streamline.  Hence 

R 


Af  =  [  -ia  (U2  +  2gA)2  (1  +  a2)2  dA  . 

'-'Re 
c 


Ii  we  set 


A  = 


R  e 
c 


■id 


this  becomes,  after  a  little  simplification 

a 

Af  =  iR  /  (U2  +  2gR  e"10)2  Jl  +  (T  e~X°  dd  .  (5.15) 

CJ0  C 

___  -|  ^ 

Here  (7  =  c(A)  =  3  .  This  method  does  not  require  that  the 

function  cr(A)  be  inverted  and  it  has  the  additional  advantage  that  it 
will  work  independently  of  whether  O’  assumes  the  value  -1  on  the  free 
streamline  or  not.  It  will  fail  if  the  free  streamline  does  not  drop 
to  a  level  where  A  =  R 
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Once  the  quantity  Af  has  been  obtained  the  depth  of  the  fluid 
far  upstream  can  be  obtained.  The  imaginary  part  of  Af  which  we  will 
call  A f  represents  the  total  flux  of  fluid  down  the  stream.  Since  the 
velocity  upstream  is  equal  to  U  ,  the  depth  h  is  given  by 


(5.16) 


In  all  the  examples  that  follow  we  will  use  dimensionless 

coordinates  obtained  by  including  a  real  positive  scale  factor  b  with 

the  dimension  of  length.  In  addition  the  dimensionless  quantity  — 

IT 

appears  quite  often  and  we  denote  this  for  simplicity  by 


F  -  £gb 
'  U2  ' 

The  position  of  the  corner  depends  only  on  c r(A)  ,  but  the  orientation 
of  the  bed  at  the  corner  and  the  depth  far  upstream  depend  on  F. 


The  depth  of  fluxd  £ar  downstream  can  be  found  from  Bernoulli 1 s 
equation  if  some  assumptions  are  made.  This  equation  in  the  form  (4.02) 
gives  the  speed  V  on  the  free  streamline  as 

V  =  [U2  +  2gAf  . 

If  we  assume  that  the  velocity  of  the  fluid  is  nearly  uniform,  then  if 
we  denote  by  d  the  thickness  of  the  fluid  perpendicular  to  the  free 
streamline  far  downstream,  we  have 


ta\i  =  dV  . 


Hence,  using  equations (5. 16;  and  (5.17) 


'  »  -JJ&  ; 


..  -:-■■■  ,  :  ,,  •'  •"  .  ■■■■ 
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d  = 


[1  +  F 


(5.18) 


It  is  assumed  here  that  A  is  large  enough  so  that  our  assumptions  are 
valid. 


Example  1  a~(A)  =  -ba/Aa  a  >  1  . 


This  choice  of  cr(A)  ensures  that  cr  — ►  -co  as  A  — ►  0  through 
positive  values.  We  have 


dA 


-ba 

a_i 

(a-l)A  1 


A 

c 


and  so  the  choice  of  a  >  1  ensures  that  this  approaches 
For  convenience  choose  c  =  oo„  Then  we  have  from  equation 


=  -iA  - 


.a 


(a-l)  X 

The  free  streamline  takes  the  form 


a-l 


.a 


x  = 


(a-l)  A 


a-l 


-  oo  as  A  0 . 
(5.02)  that 


(5.19) 


(5.20) 


This  curve  is  shown  in  figure  5*04. 


Figure  5*04 
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In  this  case  A  is  determined  by 


a 


(A)  =  -  -  =  -i  . 


That  is 


-i JS- 

A  .be'211 
c 


(5.21) 


This  has  a  negative  imaginary  part  and  thus  corresponds  to  a  value  of 
z  beneath  the  free  streamline.  In  fact,  if  we  use  equation  (5.19)  we 
have 

Ka  A_ 

—  -f  "X  ^ 

ZC  “  c  a-1  X  a 

c 


lit 


ab 

a-1 


2  (l+  l/cc) 


as  the  location  of  the  corner. 


a 


We  have  that  o'*  (A)  = 


ab 

Aa+1 


a 


cj!(AJ  = 


ab 

Aa  A 
c  c 


This  gives 
-a  -|i(1'1/a) 

r  e 


as  is 


Since  this  is  not  zero,  and  since  equation  (5.21)  shows  that  A  is  not 
a  negative  real  number,  then  the  angle  c^  of  figure  5° 03  is 
shown  in  figure  5.04.  In  addition,  if  we  use  one  of  equations  (5. 05) 
then 


c 

2n 

3 


|i  -  "2 


it 

2a 


Using  equation  (5,21)  in  the  other  of  equations (5.05)  we  bave 


v  =  tan 


-1 


-gsb  5ln  St 

LU^  +  2gb  cos  ^a  - 


=  tan 


-1 


-F  sin 


2a 


1  +  F  cos  ” 


where  we  have  made  use  of  equation  (5.I7).  These  results  and  equation 
(5.06)  give 


■  si 


j. o  :d:,  ":.o  ■.  .  oj.  -i-y 

■"  .  :  ■  \  ...  .  ,  w «•. 
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c  = 


2 1 
3 


1 

3 


r 

r  -F  sin  “  n 

Tt 

*  -  O  I-™*1 

2a 

'  2  + 

2a  "  2 

in  3t 

_ 

Ll  +  1  cos  ”  J 

2L  .  2L-  _  1  - 

2  6a  3 


F8lnfe 


1  +  F  cos 

2a 

This  relates  the  angle  c  and  F  when  a  is  given. 

The  quantity  v  of  equation  (5*05)  can  be  written 


(5.22) 


v  =  arg[l  +  F  ”] 


and  this  is  zero  when  F  is  zero.  Hence  we  choose  the  branch  of  the  inverse 
tangent  in  equation  (5.22)  to  be  the  one  which  is  zero  when  F  =  0.  If  a 
is  fixed,  then  c  can  easily  be  shown  to  be  a  monotonic  decreasing  function 
of  F.  A  sketch  of  c  versus  F  is  shown  in  Figure  5.05. 


This  shows  that  it  is  impossible  to  ever  arrange  it  so  that  c  =  0,  that  is 
that  the  bed  leaves  the  corner  horizontally.  The  bed  will  enter  the  corner 
horizontally  if  c  =  ~  .  Since  F  is  positive  this  can  be  arranged  only 


if  1  <  a  <  3- 


'-i  {C.:S,::;/  !i.v 

'  l.  :  ;i o  d:.v;  4  vi 
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We  now  examine  the  flux  of  fluid  down  the  stream.  Since  the  free 
streamline  has  a  slope  of  -1  at  some  point  for  all  a  >  1  ,  then  we  can 
use  the  integral  (5.14)  to  compute  Af.  We  have  Aa  =  -ba/cr  ,  and  so 

X  =  b(-<j)-1/a  . 

If  -cr  =  e*0  then  X  =  b  e  d0/a,  =  -  —  e  d0/a  f  and  we  get 

at)  <x 

Af  =  -  -J~2  f/  (U2-f2gb  e”10/a)'2  (cos0)2ei9/2(~)e"10/a  d0 

'o 

p*/2 

J 

0 

Two  special  cases  suggest  themselves:  the  case  where  F  is  zero  and  that 
where  it  is  very  large.  These  correspond  respectively  to  rapid  and  slow 

■austcz- 

flow  upstream. 


ib  n/~2  U 

a 


( l+F  e-10/a)4 


COS0)2 


i0(s-z) 


dO 


5*23) 


if 

Af 


F  =  0  then  we  have 


r/2(coso)i 

(X  Kj 

0 


d0 


The  change  in  the  stream  function  is,  then 

,jt/2 


Af 


-  !  (cos0)i  ^[e(§-i)]de 


J 

0 


If  on  the  other  hand  the  flow  upstream  is  very  slow  then  F 
so,  approximately 


Af 


ib  n/~2  U  n/f 

a 


pW2 

J 

0 


l.  i  §[l-3/a] 
(cos0)2  e 


d0  . 


is  large  and 


This  gives 


.  A 


imp ''‘p ■; .3 
}  ■  ■■ '  ■  i  ■■?... 


'* ;  p.'-c.'iqqp-  ,  >® 
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£aj/ 


b  n/~2  U  n/f  f*/2  ,  .^4  re  . 

=  - - -  /  (cose)2  cos [~1 1-3/a) J 

CO  ,  / 


J 

0 


de 


The  other  case  where  the  integral  for  Af  separates  easily  into  real  and 
imaginary  parts  is  when  F  =  1.  If  this  is  the  case  then 


_ie  _ie  ie 

ft  •£?  a  \~2  r  2a,/  2a, 

(1  +  F  e  )2  =  [e  (e  + 


ie 


ie 


2a)P  =  e'^  (2 


COS 


2a 


Hence 


Af  = 


ib  J2  U 


a 


J 


p«/2 

/  (2cos  @ 


2a 


cose): 


ie(i-S-) 


de 


0 


and  so 


The  values  of  these  integrals  enable  us  to  calculate  the  depth  h  of  the 
stream  upstream  and  the  thickness  d  perpendicular  to  the  free  streamline 
at  a  distance  A  below  the  origin.  These  are  given  by  the  relations 


and 


d  = 


1 

2 


The  results  are  tabulated  below  together  with  the  depth  of  the  corner  below 

the  origin  which  is  given  by  cos  r“  ,  alues  of  a  =  2  and  a  -  3 

a-1  2a 


were  used. 


'  F  = 

=  0 

F 

=  1 

F  =  cx 

Depth 

h 

d 

h 

d 

h 

d 

of  Corner 

a=2 

0.85b 

0.85b 

1.18b 

1.18b 

A  I 

0.83b  VF 

0.83b(|)"2' 

4~2  b=1.4lb 

a=3 

0.56b 

0.56b 

0.79b 

0.79b 

u# 

0.57b  n/f 

o.57t>(|F 

b=l„30b 

» 

.  1 :i-f. I'- 
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In  the  special  case  where  a  ■=  1  much  of  the  above  work  is  valid. 
Here,  if  c  is  chosen  equal  to  b,  then 


and  so 


-  f  <r(A)  dA  =  b  in(~)  , 

w 

C 

z(X)  =  -iX  +  b  in(~) 


The  free  streamline  is  given  by 

k  =  b  in(-) 

which  is  sketched  in  Figure  5«06  • 


Figure  5.06 


It  is  easily  verified  that  the  corner  corresponds  to  A  -  -ib  which  is 

c 

given  by  equation  (5 .21)  with  a  =  1,  This  in  turn  gives 


z  - 
c 


Equation  (5.22)  is  valid  when  a=  1,  The  result  is 


c  =  “  -  |  tan"L(F) 
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As  F  increases  from  zero  to  infinity,  c  decreases  from  —  to  zero. 
These  values  correspond  respectively  to  the  bed  entering  and  leaving  the 
corner  horizontally. 


The  integral  (5.23)  is  valid  when  a  =  1  and  so  the  three  cases 
F  =  0,  F  =  1  and  F  =  00  correspond  to 

oJt/2  1.  a 

=  b  v2  U  /  [cose]2  cos(-) 

0 

jt/2 


de  =  bu  |  , 


C#  =  2b U 


r 


[cos  "  cose]2  cos(*^~)  de  =  1.95  bU  > 


& 


J 


zs p 


=  V~2  bU  Vf  J  ^  [cose]2  cose  de  =  1.24  bU  v/f  . 


These  give  the  depth  h  upstream  as 


h  =  ~  ,  1.95b,  and  1.24b  n/f 


when 


F  =  0 ,  F  =  1,  and  F  =  oo 


respectively.  The  corner  is  at  a  depth  of  ~  b  =  1.5Tb  below  the  origin. 
Both  the  cases  F  =  1  and  F  =  co  have  the  bed  rising  to  the  corner. 


The  case  F  =  0  is  of  special  interest.  It  can  be  considered  as 

a  gravity  free  case  and  it  has  the  property  that  the  corner  and  the  bed  far 

upstream  are  both  at  the  same  depth.  In  fact,  the  whole  bed  between  these 

jrb 

points  is  flat  and  horizontal  at  a  depth  of  —  .  To  see  this  consider  the 
points  corresponding  to  A  =  — for  £  >  1  and  real.  This  gives  that 


z  =  -i  +  b  int^) 
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which  traces  out  the  required  horizontal  curve  as  £;  runs  from  1  to  oo. 
To  see  that  this  is  a  streamline,  we  remember  that,  if  F  =  =  0,  then, 
from  equation  (l|.ll) 


2sb 


df  =  U(l+(T  )2  dX 


We  have  cr  =  -  =r  =  -if-  and  so 


df  =  u(i  -  fy  (%  d| 


which  is  real  if  £  >  1,  Hence  the  change  in  f  with  |  >  1  is  real  and 
so  this  is  a  streamline.  The  fact  that  it  is  horizontal  is  further  born 
out  by  using  equation  (  5.0l)  with  F  =  0; 

u  =  u{2^1)2 

V(T+  i' 


w  =  U(“^7y)2  =  u(^^)2  which  is  real  for  ^  >  1, 


This  gives 

This  means  the  velocity  is  horizontal  for  all  j|  >  1 


The  thickness  d  of  the  stream  at  a  depth  A  below  the  origin 
is  given  by  equation  (5.18).  This  is  measured  at  right  angles  to  the  free 
streamline  and  so  is  nearly  horizontal  for  large  enough  A.  The  values  are 

d  =  h  =  ,  d  =  and  d  =  1.2M>(|)S 

(I+j)2 

if  F  =  0,  F  =  1 ,  and  F  =  oo 

respectively. 

Some  sketches  appear  in  figure  5«1^-  at  the  end  of  this  chapter 
which  show  how  the  angle  c  varies  with  F.  The  sketches  include  the 
three  cases  a  =  1,  2,  5.  The  depth  far  upstream  is  indicated  on  each 
diagram  except  where  it  is  infinite.  These  sketches  appear  on  page  75. 


■  . 


..  *  ■,  •'  ■'  :•  ■  ■  ■  :  .  ...  >■:.  :■  ' 


1  *  ■■  . 

.....  .• 
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Example  2. _ cr(A)  =  -cot  ft  -  b/A  0  <  3  <  fl/2 


Here  again  cr(A)  tends  to  -  oo  as  A  ->  0  as  it  should.  We 


have 


-  r  a(A)dA  =  (Acotp  +  b  in  A) 


and  so  if  we  choose  c  =  b  then  this  and  equation  (5.02)  yield 


z  =  -iA  +  (A-b)cotp  +  b  in(~) 


(5.2^) 


The  free  streamline  is  plotted  in  Figure  5.07  and  is  given  analytically  by 

x  =  (A-b)  cot(B  +  b  in(~)  .  (5.25) 


The  slope  of  the  free  streamline  tends  to  -tan|3  far  downstream,  but  the 
curve  is  to  the  right  of  the  line  x  =  (A-b)cotp  by  an  amount  b  -6n(— ). 
The  approach  to  the  x  axis  is  exponential  here  whereas  it  was  algebraic 
in  example  1. 

The  corner  here  is  determined  by  A  where 

c 


-  5^  - 


a(\)  =  -co tp  -  =  -i  . 

c 

Hence 

b  e‘lp 

sinp  * 

and  so 

A^  =  -b  sinp  =  b  sinp  e  ^  (5.26) 

as  is  consistent  with  the  result  that  Im(A  )  <  0.  This  of  course  ensures 


that  z  = 
c 

z(A^)  lies  beneath  the  free  streamline.  Since  A  (cotp-i)  =  -b , 

then 

=  A  (cotp-i)  -  b  cotp  +  b 

=  -b(l+cotp)  +  b  In  sinp  +  i  b(p-xc) 

=  -b [ l+cotp-;-incotp ]  -  i  b(^-p)  . 

This  gives  the  position  of  the  corner  in  the  third  quadrant. 


and  so 

For  this  example  cr"  (A)  -  ~2  .  Hence  we  have 

_i  /t,  \  —  A — _ -  ~2ip 

( A  1  =  2  e 

c  b  sin^p 

M  =  arg[-<r?  (A  )  ]  =  jt-23 

using  one  of  equations  (5.05).  We  again  have  the  situation  where  O'*  (A  ) 

c 


^  0  and 

A^  is  not  a  negative  real  number.  We  have  also  that 

where  F 

v  =  arg(U2+2gAc) 

=  arg(U2-2gbsinp  e*^) 

_  tan” if" _ -~F  •s-i^  -1 

*-l-F  sinp  cosp^ 

=  as  usual.  Finally  equation  (5.06)  gives 

<=  -  %  -  f[-ap-2tm-1[ir;Vi g-£.gll 

(5.27) 

.  m  .  2  ,  -1  r _ F  sjnfe.,  , 

3  3  Ll-F  sinp  cospJ  ' 

2& 

3 


Again  c  is  a  monotonic  decreasing  function  of  F  and  takes  the  value 
when  F  =  0.  Here  the  curve  is  concave  down  if  F  <  cot(3  and  concave  up 
if  F  >  cot(3  .  A  sketch  of  c  versus  F  is  found  in  Figure  5*08. 


Figure  5.08 


The  inflection  point  is  to  the  left  (right)  of 
greater  than  (less  than)  ^ 


F  = 


1 

sin(2p) 


when  p  is 


The  bed  will  leave  the  corner  horizontally  for  F  =  c  .c  ,2(3)  and 

this  is  possible  for  all  (3  in  the  range  0  <  3  <  “  .  However,  since 

<  ■“  ,  it  is  impossible  to  arrange  it  so  that  the  bed  enters  the.  corner 
3  3 

horizontally  if  0  <  3  <  —  .  If,  in  fact,  3  =  ~  we  have  the  flow  dis¬ 
cussed  in  Example  1  with  a  =  1.  This  has  the  bed  entering  horizontally 
when  F  =  0. 


In  this  example,  unless  3  >  ^  cr  will  never  be  equal  to  -1 
on  the  free  streamline,  and  so  we  cannot  use  the  integral  (5.14)  to  compute 
Af .  However  (5.15)  will  work  for  all  3.  W<  have  that 


~h 

c 


=  b  sinp  e 


-i(n-3) 


-ia 


c 


> 


=  R  e 
c 


'  •!  ,  ( : 


■■  h-/v  !  ^d:l! 


and  so 


and 
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R  =  b  sinp 


a  =  rt-B„  If  we  set 

c 


-v  -10 

A  =  R  e  =  b  sinp  e 


-10 


then 


Hence 


1  10  _  10 

/w  b  0  e  -cosB-e 

a(A)  =  -cotp  -  t  =  -cot p - — —  =  - 7~ - 

A  H  sinp  sinp 


[l+a2(A)]2  =  -^-p[sln2p-f(cos'1p+2cos|Bei0+e2ie)  ]2 


i£ 

/2  e  2 

sin(3 


[cos(3  +  cos0]2 


Finally  integral  (5. 15)  gives 

.jt~p 

0 


l£ 

1  r  2  2 

Af  =  ib  sl.n(3  /  (U2+2gbsin(3e  (cosB+cos0)^e  ^d 9 


r 


1  7t  -  P 


'10x4, 


10 


=  ib  ^2  U  J  (l+Fsinpe  xu)2(cos|3  +cos0)2  e  ^  d0  .  (  5*28) 


This  again  separates  into  real  and  imaginary  parts  easily  for  three  choices 
of  F,  If  U  Is  very  large,  then  F  is  nearly  zero  and  so,  approximately, 


rit-p  1 

Af  =  ib  v2  U  j  ( cos P+COS0 ) 2 


A 

JL  “  ip 

e  ^  d0 


This  gives 

r  P,t"P  k  q 

A^  «  b  v2  U  /  ( cos (3+COS0 ) 2  cos  -  d0  , 

kJ 

0 

If,  on  the  other  hand,  U  is  very  small,  then  F  is  large  and  we  have, 
again  approximately,  that 

X  nJt-'p  1. 

Af  e  v 2  ibU(Fsinp)2  /  (cosp+cos0)2  e  d@ 


J 


0 


Hence 
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r  1  pfl-P  1 

Alp  =  bU  v2(Fsinp)2  I  (eosp+cos 0)2  cos0  d0 

w 

0 


One  more  special  case  suggests  itself;  If  we  choose  F 
can  write: 


sin(3 


then  we 


ie 


(l+Fsinpe  i0)2  =  e  H  (2cos  ^)2 


0< 


Hence  equation  (5.26)  becomes 

i  Jt-(3 
0 

Here  then, 

i  Jt-p 


Pit*  P  a  X  ~  h 

Af  =  2ibU  I  [cos  —  (cosp+cos0) ]2  e  d0 


p7iMp  Q  A  Z2LQ 

Ap  =  2bU  I  [cos  —  (cosp+cos0)  ]2  cos  ^ 

'  0 


d0 


when  F  = 


sing  * 


As  in  Example  1,  these  integrals  were  evaluated  for  special 
values  of  p.  The  values  chosen  were  p  =  g  ,  p  =  ^  and  p  =  ~  . 

These  results  were  used  to  compute  the  depth  h  far  upstream  and  the  thick¬ 
ness  d  at  great  distances  downstream.  The  corner  in  this  case  is  at  a 
depth  ( it — P ) b  below  the  origin.  All  these  results  appear  in  tabulated 
form  below. 


F  =  0 

F  =  l/sinp 

F  =  00 

Depth 

of  Corner 

h 

d 

h 

d 

h 

d 

p  ‘  6 

2.93b 

2.93b 

’  3.06b 

3.06b 

[uf]i 

1.05b  n/f 

1.05b(|)2 

'  |pb=2.62b 

P  ~k 

2.69b 

2.69b 

2.92b 

2.92b 

r  1  12 

lu—P 

1.30b  Jf 

1.30b(|)2 

i~h-2o 36b 

P  =  3 

2.56b 

2.36b 

2.68b 

2.68b 

2A  X 

‘1+7P2 

1.42b  Sf 

1.42b{^)i 

•“b  =2. 09b 
3 

•  -V 
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Sketches  of  the  case  (3  =  jj-  appear  in  Figure  5«!5*  They  give 
an  indication  of  how  the  angle  c  varies  with  changes  in  F„  Where  h 
is  finite,  it  is  indicated  on  the  diagram.  These  sketches  appear  on  page  76, 


2 

Example  3.  cr(A)  =  -cot  (3[1  +  ^75].  0  <  |3  <  x/2. 

In  this  case  we  have  that 

-  r  a(A)  dA  =  cot(3  [A  -  ^r-]^  =  cot(3  [A  -  ^r-] 

J  A  c 

c 

if  we  choose  c  =  b.  Hence,  using  equation  (5.02),  we  have 

z  =  -iA  +  cotp[A  -  “]  .  (5.29) 

The  free  streamline  is  given  by 

x  =  co t (3  [A  -  “]  ,  (5.30) 

and  is  sketched  in  Figure  5«09.  If  A  ^ 

“b2  cotB 

x  *  A 

while  if  A>  >  0,  we  have 


x  -v  A  cot  3 


:i:  r:i 


~-;d!  '  ;  a  :;  • 
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Hence  the  curve  approaches  the  x-axis  and  the  line  x  =  A  cotp  algebraic¬ 
ally.  The  slope  of  the  free,  streamline  at  the  point  (x,y)  =  (0,-b)  is 

tan(3  . 


We  have  A  given  by 


o-(Ac)  =  -cotp  [1+^]  =  ~i  . 

This  gives  v— 5  =  i  tanp-l  =  ~~r  e  .  Hence 

cosp 

A2  =  -b2  cosp  e*^  , 


(5.51) 


mm 

Choosing  -1  =  e  so  as  to  make  Im(A  )  <  0,  this  giv- 

r - 

A  =  b  vcosp  e 


(502) 


Taken  with  equation  (5.29)  this  gives,  after  some  algebra,  that 

zc  =  Mcot  P"1)  "  b2cotP  ) 

c 

_  2b  n/cosP  e'1^2+l^ 
sin  p 

is  the  location  of  the  corner. 


We  now  calculate,  the  angle  c,  We  have  that 


and  so 


<x(A)  =  -cotp  [1  + 

8  / *\  ^  2b2  cotp 

a  (A}  *— ^3  • 


Hence  using  equations  ( 5 . 31 )  a*id  (5«32)  and  some  algebra  we  get 

,  ^  2b2  cotp 

<V  -  X"  'V>  =  ~~ 

c  c  sinp 

Hence  the  first  of  equations  (5.05)  gives 


-2 


i(0&) 

M2  2 ' 


-s/cosp  b 


N 


....  Jr.  . ■.>./  -  . .  v 
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H  =  arg(-a'(Ac))  =  |  -  ^ 
The  second  of  equations  (5*05)  gives 


v  =  arg(U2  +  2gA  ) 


arg(u2  +  2gb  \fcos^  e 


=  tan 


-1 


■F  n/cos(3  cos 


-1  +  F  \/cos(3  sin  ^ 


Ppb 

where  F  =  .  These,  taken  with  equation  (5.0 6)  gives 


Jt 

1 

it 

2tan  1 

“F  vcos(3  cos  g  - 

- 

u  "  3  ’ 

3 

2 

2 

-1  +  F  n/cosB  sin  ^  - 

- 

I  JJ 

6  2  3 


tan 


F  n/cosP  cos  2  ‘ 
1  +  F  n/cosP  sin 

2 


(5.33) 


This  gives  c  as  a  decreasing  function  of  F.  The  curve  is  sketched  in 
Figure  5.10. 


The  range  includes  c  =  0  if 
the  corner  horizontally  only 


Tt 

6  <  —  and  so  the  bed  can  be  made  to  leave 

“  5 

if  this  condition  is  satisfied.  The  value 


. 


-1  '  -  '  *  ,  v, 

.0.  .?  ;'X«4:i7 


1  5-'  v  Jfc  t/  •:  '  C  ",  '•  :  rr  i.  ■  •= 
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T(  jC 

c  =  —  can  occur  only  if  p  >  —  .  Hence  this  is  the  condition  which  must 

3  ”3 

be  met  in  order  that  we  can  choose  F  so  that  the  bed  enters  the  corner 
horizontally. 

Here,  as  in  Example  2,  cr  need  not  necessarily  equal  -1  at 
some  point  of  the  free  streamline.  However  we  can  again  set  up  an  Integra' 
for  Af  in  the  A-piane  using  equation  (5»15)*  We  have  that 


/ - -  Mp  p/ 

A  =  b  vcosp  e 


and  so  R  =  b  \]  cosp  and  cc  =  ^  "  p  •  If  we  a  gain  write 

0  C  di  cz. 


'i6 


then 


A  =  R  e  ^  =  b  n/cosB  e 
c 


cr2  =  cot2p  [1  +  +  --] 


■19 


**  T  AV 
[cot2p  +  2c.osp  e2l9+  e^it'  J  > 


using  the  equation  <r(A)  =  -cotp[l  +  ^2  J .  It  follows  that 

A 

/T“”2  1  M  o  *  210  4ie,4 
vl+c  =  TT^T  [1  +  2cosp  e  +  e  j2 

S  Xup 


v/2  e1®  ,4 

=  ---^7™  [cosB  +  cos 29 3 _ 

smp 


Hence  equation  (5.15)  becomes 

Af  =  ib  Vcosp  I  2  ^  (U2+2gb  n/c asp  e 
0 


eL& 
s  in  3 


JL  „  i  g 

(cosp-:-cos20)2e  ~  d@ 


LJk±Ji£21&Jl  j  2  2  (1+F  %/4osjJ  e"lf9)2(cosp+cov2G)2  d0 


sin  p 


J 


0 


(5.3*0 


where  F  - 

U 


- 


i  '  . 1  '  - 
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If  U  is  very  large,  then  F  is  nearly  zero  and  we  can  write, 
approximately,  that 

Af  =  f  2'2  (COS|3+cos2e)i  d 9  . 

sin0  J  ' 

0 

Hence  the  change  in  the  value  of  the  stream  function  from  the  bed  to  the 
free  streamline  is 

=  b~U's^p°S^  J'Z  2  (cosp+cos20)2  d0  . 

'  0 

In  the  case  where  U  is  very  small,  F  is  large  and 


Af 


It 

_  iU  J~2  b  \/f(  cos(3  rl-f 

sinp  j 

0 


W 

A  ~  O 

(cosp+cos2@)2  e  d Q 


2L.i 


Hence; 


f  (cosp+cos2@)2cos| 


As  in  Example  2  there  is  one  other  value  of  F  which  allows  us  to  split 
the  integral  in  equation  (5*37)  into  real  and  imaginary  parts.  If 


F  = 


\! cos  3 


then 


-i- 

( 1+F  slcosfi  ei0)2  =  e  ^  (2  cos  —)2  , 


and  so 


jt_J3  _±9 

Af  =  ~ibsinpS~^~^  2  e  ^  tcos  §(cosP+cos2@)  ° 

0 

_ 

Hence  A tf/ =  ^"sin^5^  2  [cos  ^(cos0+cos2@) ]2  cos  ^  d0  , 

0 


when  F  = 


n/cos(3 


I 


■’  .  :  ;  ,  ;t. 

*  7  .  .  ‘  i  . 

.1  1  cr -.-rS'-lv:  Sfj.r*  J 


:  /■ '  ' 


:  -V-M  '  '  ■  ..  -  ■  :  ■  V"  <  .  i,  ;  r.  ' 


iTTt  *u  V  .  . 


M 


;■  •  vs 


I  tv  ;!■;? 
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As  in  the  above,  these  integrals  for  A  if/  have  been  evaluated 
for  p  =  —  ,  ^  .  The  depth  h  upstream  and  the  thickness  d  down¬ 

stream  perpendicular  to  the  free  streamline  are  tabulated  below,  together 
with  the  depth  below  the  origin  of  the  corner.  The  latter  is  given  by 

2b  VcosB 


sinp 


cos 


F  =  0 

F  =  l/</cosp 

F  =  00 

Depth 

of  Corner 

h 

d 

h 

d 

h 

d 

Jt 

P  =  6 

3.31b 

3.31b 

U.  54b 

4.54b 

r,  -fz-K 

-  V~3  b 

~T 

2 

3.04b  n/f 

3.04b(|p 

3.60b 

0  =  ~ 

P  4 

1.88b 

1.88b 

3.06b 

-.■,3.,0.$b 

r.  Vz^i 

L  b  J 

r 

2 

1.65b  n/f 

1.65b(|F 

2.12b 

n  11 
p  =  5 

1.10b 

1.10b 

l.52b 

1.52b 

-  SzXl 

+  b 

1 

2 

0.89b 

0.89b(^F 

1.41b 

Example  4  <x(A)  =  — •  0  <  (3  <  jt/2. 

Here  again  (j(A)  -*■  -  oo  as  A  ->-0.  In  this  case  we  have 

-  f  or(A)dA  =  cotp  J'  —  dA  =  b  cotp[!n( e^^-l) 

c  '  c  e  "1 

If  we  choose  c  =  b  /n2,  then  in(eC/^>-l)  =  0  and  so 

z(A)  =  -iA  +  b  cotp  .^[e^/k-l]  »  (5.35) 

The  equation  of  the  free  streamline  is 

x  =  b  cotp  /nfe^^-l] 


(5.36) 


■ 


•*s  7i"-;  '  "'  f  •  '  v- -Vj  ;:v, 


O  ii L.'  !;! j; 


a--. 


-  6k 


If  A  >  >  0  then  x  2;  b  cot (3  ( )  =  A  cot  (3  .  If  ,  on  the  other  band,  A  *\j  0, 

a  .*;.fnP  x 

then  x  ^  b  cot  3  in  (— ),  or  A  <v  b  e  b  ,  This  is  an  exponential  func¬ 

tion  for  x  <  <  0  and  so  the  approach  to  the  negative  x-axis  is  exponential. 
In  addition,  if  we  write 

x  =  b  cotp  in[e  ¥>  (l-e“M>)] 

=  A  cotp  +  b  cotp  in( 1-e  , 

then  for  large  A,  x— Acotp  ^  ~b  cotp  e  which  shows  that  the  free 

streamline  approaches  the  line  x  =  A  cotp  exponentially.  The  slope  of 
the  free  streamline  is  -■§■  tanp  at  the  point  (x,y)  =  (0,-b  in  2)„  The 
curve  is  sketched  in  Figure  5.11. 


Figure  5 .11 


In  this  example,  A  is  given  by 


<V  ■  = 


This  gives  -i  cot  p  =  1  -  e 


-\> 


and  so 


ie 

e  =  1  +  i  cot  p  = 


-ip 


sin  p 


(5  37) 


v.'iXS 


■  1  ■-  '  -tj 


;  :'7'  I  ; 


jsa. 


3  •'  I_T  ■, 


& 


Finally 
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-  in  sinp  -  i(|  -  3)  . 


(5.38) 


To  compute  the  position  of  the  corner,  z  ,  we  first  calculate  e  -  1, 

c 

Equation  (5*37)  gives 

\A> 


and  so 


e  C  -  1  =  -i  cotp  e  C 


A  Ay 

in[e  -  l]  =  in  cotp  +  -  i  ” 


=  in  cosp  -  i(jt  -  (3) 


using  equation  (5*38).  Finally,  then,  equation  (5.35)  yields 


z 

•—  =  -i [in  sinjB  -  i(~  -  p)]  +  cotp[in  cosp  -i(?t  -  p)] 
=  "[(f  ~  0)  +  cotP  secp3"i[in  sinp+cotp(jt  -p)]. 

(5.39) 


Since  in  this  case  er(A)  =  — ,  then 


1-e 


o-’(A)  = 


cotp(-e  A//^)(  £) 

[l-e-AA]2 


-  cotp _ 

b[l-e'AA][i_eAA] 


We  have  from  equation  (5*37)  that 


and 
and  so 


CT'(X  )  = 


V, 

[1-e  ]  =  -i  cotp 

ACA 

[1-e  ]  =  i  cotp  e 

_ -cotp _ 

—f— 

b [-icotp] [icotpe  C  ] 


V. 


cf  b 


Kf-3) 


-  tan  P  e 

b  sin  P 


using  equation  (5.38).  Hence 


'>  a.J  ff'3  Jfl  ;  a 


Jbraa 
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In  addition 


p  =  argC-CT'  (Ac))  =  |  -  P  . 


v  =  arg[U2  +  2gAc] 

=  arg[U2  +  2gb ( in  sinp  -  i(|-p))] 


tan 


-i  r  -f  (5  -  3) 


F  in  siup 


where,  again,  F 
obtain: 


=  .  Using  equation  (5.0 6)  with  these  results  we 


Jt  1 

jt  „  _  ^  -1 

r  -F  (|  -  p)  -1 

-| 

C  '  3  "  3 

-  -  p  -  2  tan 

1  +  F  in  sinp 

jt  £  2  . 

=  6  +  3  '  3  ta" 


-i  r  F(l '  p) 

1  +  F  in  sinp 


As  in  previous  examples  this  is  a  decreasing  function  of  F  if  p  is 
fixed.  The  angle  c  vanishes  when 


F  =  [(|-p)  cot(^+|)  +  in  esc  p] 


Since  the  right  hand  side  is  positive  for  all  values  of  p  in  the  range 
0  <  p  <  ~  ,  then  this  shows  that  for  each  p  there  exists  an  F  which 
makes  c  =  0,  i.e.,  which  makes  the  bed  leave  the  corner  horizontally. 


When  F  =  0,  c  =;+f  .  When  F 

o  5 


£  .  E  _  2 
3  6  3 


tan 


i  - g 

in  cscp 


_  £  ^  E  £ 

3  6  3 


jt  “  tan  Y  2  ^ 


in  cscp/ - 


- 1  •  I +  f can' 


E  . 


P 


in  cscp 


where  the  inverse  tangent  in  the  last  expression  has  its  principal  value. 


-  !■ 


c  .;y.  •  ‘  ", 
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then  it  is  impossible  that  c  = 


Since,  for  given  P,  c  <  ^  +  ■§ 

o  3  j 

unless  p  =  ~  .  That  is  to  say,  for  any  p  <  —  ,  it  is  impossible  to 

Ti 

choose  F  so  that  the  bed  enters  the  corner  horizontally.  If  p  =  —  , 


then  F  must  be  zero  in  order  that 


n 

c  =  1 


Example  5  cr(A)  =  -sec2  (—  - 


This  example,  unlike  the  others,  has  the  property  that  the  free 

p  A 

streamline  is  bounded  in  the  y-direction.  We  have  -  /  cr(A)dA  = 


J 


f  sec2(g  -  ~)dA  =  b  tan(j^  -  ~)  if  c  is  chosen  equal  to  ™ 


This  give 


z  =  -  iA  +  b  tan(~  •  . 

b  d 


(5.40) 


The  free  streamline  is  sketched  in  Figure  5.12. 


The  equation  of  the  free  streamline  is 

x  =  b  tan(£  -  |)  . 

The  corner  of  the  flow  is  given  by  A^  where 


(5.41) 


, 


■;  ‘  '  '  .--r'  ?;C 
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A 

r_i 
k  b 

A 


a(Ac)  „  -  sec2{~  -  f)  *  -  1  , 


-hi  2 1 


cos2(f  -  f)  -  e 


C09<r  - 1>  -  ±  (A  -  ^ 


If  we  write  (—  -  ~)  =  p  +  iq,  then  if  we  choose  cos( 


c  n 


J_  A. 
s/2  '  42 


we  have 


_L  __i 

cos  p  cosh  q  -  i  sin  p  sinh  q  =  ^ 


Hence 


This  gives 


i 

cos  p  cosh  q  =  j-g  =  sinp  sinh  q  . 

“  =  cos2?  cosh2q  =  ( l-sin2p)(  l+sinh2c 


| 

If  we  use  the  fact  that  sin2p  sinh2q  =  —  ,  this  gives 


and  so 


sin2p  =  sinh2q  , 

k  4  1 

sin  p  =  sinh  q  =  “ 


These  yield  P  =  ±  0.999>  ft  ±  0.999 

q  =  +  O.765  . 

Since  we  want  Im( A^)  <  0,  then  q  =-0.765.  Using  this  results,  the  two 
conditions  cos  p  coshq  >  0  and  sin  p  sinh  q  >  0  imply  that  cos  p  >  0 
and  sin  p  <  0.  This  means  that  p  lies  in  the  fourth  quadrant  and  so 
p  =  -0.999*  Finally 

A 

("jf  -  |)  =  -0.999  "  0.765  i  , 


A 

~  =  0.572  -  O.765  i  . 

b 


In  the  other  case  where 


A  «.  ^+-^- 

cos(T  ’  a>  =  Sz  /2 


cos  p  cosh  q  -■  i  sin  p  sinh 


q  =  nT2  V2 


we  obtain 
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2l 

Hence;  cos  p  cosh  q  =  =  sin  p  sinh  q 

This  leads  to  the  same  possibilities  for  p  and  q  and  again  Xm(A^)  <  0 
requires  q  =-0.765.  However  in  this  case  the  relations  sin  p  sinh  q  <  0 
and  cos  p  cosh  q  <  0  give  that  p  lies  in  the  second  quadrant  and  so 
p  =  it  -  0.999.  Thus: 

A 

(■f  -  f)  =  Jt  -  0.999  -  0.765  i, 

A 

=  *  +  0.572  -  0.765  i  . 


We  can  now  obtain  the  position  z^  of  the  corner.  We  can  write 


that 


tan(”  -  |)  =  ±[sec2(-~  -  |)  -l]2  =  +  (i-l)2  =  +  V2  e  ^ 


That  is 


tan(~  -  |)  =  +  (0.455  +  1.0991)  . 


A 


Since  Im(A  )  <  0,  then  Imj^ta^T™  -  ~)  ]  <  0  and  so 


tan(~  -  |)  =  -0.455  "  1.0991 


If  we  use  this  and  the  two  values  of  __c  in  equation  (5.40)  we  get 

b 

zc  =  (-1,220  -  I.67I  i)b  , 
and  zc  =  (-1.220  -  (I.67I  +  *)i)b  . 


We  must  choose  the  highest  value  and  so 


zc  =  -b  (1.220  +  I.67I  i)” 


This  corresponds  to 


and' 


Ac  =  b(0. 572  -  0.765  i) 

5*1 
V2  e  4 


A 

_C 

b  ~~2‘ 


tan(-^-|)  =  - 


' 


q 


■■■  38  5a:ib:0  .  i 


-  TO  - 


To  calculate  the  angle  c  in  terms  of  F  we  proceed  as  follows 
2  A  tx 

We  have  cr(A)  =  -sect“(—  -  — )  .  Hence 


i  n  \  2  2/A  it \  /A  tc \ 

a  (A)  ,  -  -  sec  (f  -  -)  tan(j  -  -) 


and  so  using  results  stated  above 


-a'(AJ  =  l  (i)Va 


)TT1 

8 


.5A  ¥ 


This  gives 


In  addition 


p  =  arg(-cr,(Ac))  =  ^ 


v  =  arg[U  +  2gAc] 


=  arg[U2  +  2gb( re -0.999  -  0.765i)l 


"  tan  I  1+2,143?  ^ 


where  F  = 

U2 


Finally,  then,  using  equation  (5.06), 


c  =  3  ‘  3  ^  '  2tan 


-1 


-0, 


1+2.143F 


2651  1  1 

"i  :zT? 


.  £  tan-i  f  0-T6?F  1 
24  3  1+2.143F1 

Here,  c  =  when  F  =  0,  and  c  -»  -O.O98  =  -(5°37?)  when  F  -*■  00 
It  is  clearly  impossible  to  choose  F  so  that  c  =  —  ,  i . e . ,  so  that 
the  bed  enters  the  corner  horizontally. 


Example  6  <t(A)  =  -  y 


This  example  differs  from  the  above  in  that  A  =  0  corresponds 
to  a  finite  point  and  that  the  flow  at  that  point  is  not  uniform.  We  have 
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if  we  choose  c  =  0.  Hence 

z  =  -iA  +  2  >/bA 


The  free  streamline  is 

x  =  2  \/bA 

which  is  a  portion  of  a  parabola  opening  down  in  the  fourth  quadrant*  A 
sketch  is  given  in  Figure  5.13* 


The  corner  occurs  at  A  where 

c 


c 


This  gives  ■n/Ac  =  -i  \Tb  , 

A  =  -  b 
c 

This  is  an  example  of  a  flow  where  it  is  possible  that  U2  +  2gA  =  0. 

G 

The  position  of  the  corner  is 


'  f  ,«0  at  'J  0  SCtOff  3  SV  3  l 

■  .  '  m<  -:i.  .  :■■■■  ■ 


^  •  ... .•</ hh  9>  ,»p>  a..l  ■ 
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=  -i(-b)  +  2(-ib)  =  -ib  . 


Also,  cr'(A)  =  ~  X  3/^  ,  and 


so 


pi  =  arg( -cr *  (Ac))  =  arg 


r  -  Vb  “ 

.2  V 


ZL 

2 


In  addition 


v  =  arg(U2  +  2gAc)  =  arg(U2  -  2gb)  =  arg(l-F) 


In  this  example,  three  cases  present  themselves.  We  shall  want  to  examine 
the  velocity  on  the  negative  y-axis  between  the  origin  and  =  -ib,  and 
to  this  end  we  make  the  change  of  variably 


This  gives 


A 

b 


=  “X 


z  =  ib|(?-2) 


which  takes  on  the  required  values  when  0  <  |  <  1.  Moreover,  if  we  use 
equation  (5.01),  then  we  can  write 


w 


(U2  +  2gA)i 


fr(A)-l\  i 
ff(A)+i ) 


u(  i+Fr) 


-i 


b 

A 


+i 


i 


Case  (a).  F  <  1. 


In  this  case  v  =  arg(l-F)  =  0  and  so 


-  %  -  5("-2v)  3  ! 


Moreover,  u>  is  real  for  all  j*  in  the  range  0  <  I-  <  1,  with,  of  course, 
a  singularity  at  |  =  1.  Hence  the  velocity  is  horizontal  on  the  negative 
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y-axis  between  the  origin  and  the  point  z =  -ib.  We  could  thus  continue 
the  flow  into  the  third  quadrant  by  reflection  in  the  y-axis  and  reversal 
of  the  direction  of  flow.  This  would  create  a  symmetrical  flow  over  a  sharp 
corner  of  interior  angle  —  with  a  parabolic  free  streamline. 

3 

Case  (b).  F  =  1. 


This  is  the  case  where  U2  +  2gA  ^  =  0.  Equation  (5*09)  gives 
2  >  and  equation  (5.12)  gives,  with  7=1, 


c  = 


_  z±i 


±JL 


tc  =  0  . 


In  this  case 

U>  =  u(i-|2)i(^|)i  =  U(i-t)  , 

and  so  the  velocity  is  again  horizontal  for  0  <  ^  <  1.  However,  in  this 
case  a  stagnation  point  occurs  at  z^  =  -ib.  Again  we  can  continue  the 
flow  into  the  third  quadrant  by  reflection. 


Case  (c),  F  >  1. 


Here  the  velocity  is  horizontal  if  o  <  i  <  and  is  vertical 

if  P  >  ~r~  *  The  point  £  =  hr  is  a  stagnation  point  and  so  the  flow 

s  vF  's/F 

cannot  be  analytically  continued  past  it.  This  singularity  corresponds  to 


z  =  y  (2nTf  -  1)  . 

It  corresponds  to  a  A  value  given  by 

fK  =  -i  _JL  . 

VTg 


Hence  U2  +  2gA  =  0  for  this  singularity  and  the  situation  is  dealt  with 
in  the  material  leading  to  equation  (5.12).  The  situation  here  has 


' 

■ 

...  ■  .(  Cf  !  i--'  <J 

- .% 

:  f  y.i'  "x  i.;".  fc.r 

, 

■  '  -  '  ■  b  . :.  i:  x  ■  .  v  x/: 


. . . ,  ■  ;0.  vix-": 


r'“  1  '  .  '  <  '  Ij. 
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which  gives 


UJ 

4T& 


Mj.  =  arg(-(j,(Ac))  =  |  . 


Since  7=1,  then  equation  (5.12)  gives 

c  =  0 

as  in  Case  (b)  above.  Again  the  flow  can  be  reflected  in  the  y-axi 


Alpha  =  3 
Figure  5„lb 


a/Aa 


Example  1  a  ~ 


M/f 
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APPENDIX  ONE 


A  THEOREM  ON  HARMONIC  FUNCTIONS 


Theorem:  Given  two  functions  u(x,y)  and  v(x,y)  which  are  conjugate 


harmonics  in  the  open  region  y  <  0.  If  u(x,y),  v(x,y)  and 

)v 


du  ,  ,  Ov 

^  exist  and  are  continuous  in  the  closed  region  y  <  0,  then  ^  exists 

Su 

and  equals  in  the  closed  region  y  <  0. 

OX 


Proof;  The  theorem  is  true  in  the  open  region  since  u  and  v  are- 

conjugate  harmonics  there.  Let  (xq,0)  be  on  the  boundary  and 

let  k  <  0  be  a  real  quantity.  Then  v(x,y)  exists  and  is  continuous  in 

c)v  du 

the  region  x  =  x^,  k  <  y  <  0,  and  ^  exists  and  equals  ^  in  the 
open  region  x  =  x^,  lc  <  y  <  0.  Hence  on  the  line  x  =  x^  the  mean  value 


theorem  gives 


v<Vk)  •  v(x0’0)  - k 


<xo>?) 


where  k  <  |  <  0.  This  is  valid  for  all  k  <  0  and  so 


lim 
k  — >  O’ 


The  limit  on  the  left  is  v~ 

dv 


f  v(x0,k)-v(x0,0)-i 


=  lim 


dv 

57 


<V 


^V 


(xQJ0) 


k 

if  it  exists.  We  know  that  inside. 


the  fluid  ^  ^  by  hypothesis.  Hence,  since  (|  -*•  0  as  k  -+0,  we 


have 


lim 

k  ->  0  L 


r  v(x0,k)-v(xQ,0)  -| 


=  lim 

J  i  -*-0  L 


Bu 


(x0,4)' 


Ou 


<V 


since  is  continuous  up  to  y  =  0  by  hypothesis.  Hence  the  limit  on 

ox 

u  in  the  region  k  <  y  <  0  . 


the  left  exists  and  equals 

ox 


<v0) 


»  ' '  ..  V‘..t  '■  ;■  ;  .  , 

fis<ja  v.ci  «|  xiitord 

2  *sis  r-  ■■■•■;•  bIvi  ■ 
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APPENDIX  TWO 

THEOREMS  ON  EXISTENCE,  UNIQUENESS  AND  CONTINUITY 
OF  SOLUTIONS  TO  ORDINARY  DIFFERENTIAL  EQUATIONS 


Lemma;  If  g(x,fi)  is  continuous  in  the  closed  bounded  region  R  given 
by  lx_x0i  1  M  and  |4-p0|  —  N>  t1lien 

f(x,n)  = 


is  uniformly  continuous  in  R„ 


f 


g(  t ,  (i)dt 


Proof:  Since  g(x,p)  is  continuous  in  the  bounded  closed  region  R, 

then  it  is  bounded  and  uniformly  continuous  there,,  The  first  of 
these  assertions  means  that  |g(x,p)|  <  k  in  R;  the  second  means  that, 
given  e  >  0,  we  can  find  a  5  >  0  such  that 


|g(x,|i)  -  g(x# ,p'  )  |  <  € 

for  any  two  points  (x,p)  and  (x',|i')  in  R  satisfying  the  conditions 

|x-x*|  <  5  and  |  p.- p. J  j  <  5 


We  can  write 

f(x,p)-f(x'  ,p'  ) 


r 

j 


[g(t,p)-g(t,|i,)]dt 


x 


0 


Px' 

J  g(t,p')dt  • 


If  ||i-|is|  <  &  this  means 

|f(x,p)-f(x\n')|  <  e  |  x-xQ  |  +  |x'-x|k 


<  e  M  +  5k 


v.rnr. 


;  ;  ■  ,'5t' *  '  imj  a 


'■  ■■■  ■  :  x  V  ; 


"  ^  ‘  7X  . : . 
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Since  this  can  be  made  arbitrarily  small,  the  result  is  proved. 

Using  this  we  can  prove  the  following  two  theorems.  The  first 
one  is  a  generalization  of  section  5 .6  and  5.7  of  [5]. 


Theorem  1;  Given  functions  f(x,y,p)  and  y^(p)  which  are  continuous 

in  their  arguments  in  the  region  R  of  x-y-p  space  given 


by 


|x-x  |  <  h,  |y-y  |  £  Mh  and  |n“kJ  <  k  .  (A.Ol) 


Suppose  that,  in  the  region  R,  f(x,y,p)  satisfies 


and 


f (x,y,p)  |  <  M 

f(x,y,p)-f(x,y*,p) |  <  N|y-y*| 


(A. 02) 
(A. 03) 


Then  the  solution,  y(x,p)  of  the  differential  equation 


S  -  £(x>y-n) 


(A. Ok) 


which  satisfies  the  boundary  condition 


y(xQ,p)  =  y0(n) 


(A. 05) 


exists,  is  unique,  and  is  uniformly  continuous  in  its  arguments  in  the 

^f 

region  R  .  Condition  (A.03)  is  met  if  ^  is  bounded  in  R  . 


Proof: 

(l)  Existence:  If  y(x. p)  satisfies  equation  (A. Ok)  then  a  quadrature 


gives 

rx 

y(x.p)  =  yQ(p)  +  '  f [x,y(x,p),p]dx 


(A. 06) 


as  a  function  which  satisfies  both  (A. Ok)  and  (A. 05)  .  We  now  define  a 
set  of  "successive  approximations"  to  y(x,p)  by 


;T'jj  •  f  ■■■■■ 


■  n 


.(W  st  wo; 


:'■>  ""  /  >  ■:  K  ;■■■ 
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yQ(x»n)  =  yQ(n) 


(A. 07) 


rx 

yn(x,Ji)  =  yQ(n)  +  j  f  [x,yn_1(x,(a),|i]dx 


(A„08) 


We  will  show  that  the  sequence  of  functions  yn(x,p.)  defined  in  this  way 
converges  uniformly  to  a  unique  function  y(x,|r)  satisfying  (A.OJ+)  and 
(A. 05).  We  use  the  notation  y^  =  y^(x,|a). 


We  first  show  that,  for  all  n. 


|yn  -  y0i  <  ®  • 


(A. 09) 


This  is  obvious  for  n  =  0  by  equation  (A.07)  and  holds  for  n  =  1  by 

conditions  (A.Ol)  and  (A, 02),  The  proof  proceeds  by  induction  on  n.  If 

equation  (A„09)  is  true  for  n  -  1,  then  equation  (A.Ol)  gives  that  y  . 

is  contained  in  R  and  so  jf(x,y  <  M  by  (A, 02).  Hence  ( A  „  C-8  ) 

n  *  i. 

gives 


|y„  -  y0i  =  J  f(x,yn.1,4)dx 


<  M  x 


This,  with  (A.Ol),  establishes  (A, 


We  now  show,  again  by  induction,  that 

.-i  lx“xoln 


I y  -  y  , j  <  MNn' 
,yn  n-11 


n 


(A. 10) 


If  n  =  1  this  follows  from  (A.Ol)  and  (A„09)°  Assume  it  true  for  n  =  p 
Now,  (A. 09)  shows  that  y^  and  y^  1  are  contained  in  R  and  so  (A.Oj) 
and  the  induction  hypothesis  imply 


f(x»yp»n)-f(*,yp.1»n)l  <  Nlyp~yp.il  <m 


x-x. 

p  ' _ o1 
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Using  this  with  equation  (A.08)  we  have 


which  establishes  (A.lO). 


If  we  now  use  (A„Cl)  with  (A0lo),  we  obtain 
.  i  „  M  (Nh)n 

IVVll  <ii^: 


which  implies  that  the  series 

00 


yo  + 


[Vyn-l! 


(A. 11) 


n=l 


converges  uniformly  to  some  function  y  =  y(xsp)  in  the  region  R,  It 
remains  to  show  that  y(x,n)  is  a  solution  of  the  system  (A„04)  and  (A„05)o 


Since  the  sum  of  p  terms  of  the  series  (A. 11)  is  y 


P-1 


then 


we  have 


lim  [y  , ]  =  y(x,n)  .  (A. 12) 

p  -4.00  P-i 

This  and  (A. 09)  ensure  that  |  y ( x ,  |_i)  ~yQ  I  <  Mh  and  so  y(x,(j.)  is  in  the 
region  R.  If  we  write  y  =  y(x,|ji)  then  equation  (A. 03)  gives 


I  f (  x.y.n)  '  f(x,yn_1,|i)|  <N|y-yn_1| 


Using  (A.08)  we  can  write 


■:  '-VS  ;  ;  .  v 


-  G3  - 


y-y 


-  r 


0  J 


f(x5y,^)dx  =  y-[y 


rx 

rj  £(**yn-i^)dx] 


r 


J 


f(x,y,p)dx 


0 


0 


=  y-y  - 

n 


f  [f(x,y,4)-f(x,yn 


1 


,n)]dx 


x. 


Finally,  then 


rx 

|y-y0"  j  f(x,y,|a)dx|  <  |y-yn|  +  N|y-yn_1|h 


using  (A. 01).  Equation  (A. 12)  implies  that  the  right  side  can  be  made 
arbitrarily  small  if  n  is  large  enough.  Hence  the  left  side  must  be 
zero  and  so  the  function  y  =  y(x,p)  satisfies 

oX 

y(x,p)  =  yQ(p)  +  J  f  [x,y(x,p),ii]dx  . 


It  is  easy  to  verify  that  this  satisfies  (A.OU)  and  (A, 05)  and  so  the 
existence  of  a  solution  to  this  system  has  been  proved. 


(2)  Uniqueness :  To  prove  uniqueness,  we  suppose  that  y(x,p)  and 
both  satisfy  the  equations  (A.04)  and  (A. 05)  where  the  conditions  (A.02) 
and  (A. 03)  are  in  effect.  These  functions  both  satisfy  equation  (A.06)  and 
so  we  can  write 

y-1  =  f  [ftx.y.nj-ffx.ri.iijldx  . 

'xo 

In  addition,  (A.06)  implies  that  the  points  (x,y,p)  and  (x,q,|j.)  are  con¬ 
tained  in  R,  and  so  (A.03)  yields 


|f(x,y,n)-f(x,n,|a)  |  <  N|y-ri|  . 


V 

'  ,  5ii rsaaiits, 
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These  results  give 


|y-n|  <  N  f  |  y-T)  |  dx  . 

d 


(A. 12) 


x 


0 


If  the  function  e 


-Nx 


«  x 

Jj 


y-rjjdx  is  differentiated  with  respect  to  x, 


then  equation  (A. 12)  shows  that  it  does  not  increase  with  x.  Since  the 
function  is  zero  when  x  =  x^  and  is  non-negative  then  it  must  be  ident¬ 
ically  zero.  This  implies  that  y(x,p)  =  n(x,|i)  in  R,  which  proves 
uniqueness . 


(3)  Uniform  Continuity;  Since  y^(p)  is  continuous  in  the  closed 

bounded  region  R,  it  is  uniformly  continuous  there.  Moreover,  since 

f  is  continuous,  then  f(x,j^(p)  ,  p)  is  a  continuous  function  of  x  and 

p  in  R  and  so  the  lemma  and  (A.08)  imply  that  y.^(x,p)  is  uniformly 

continuous  in  R.  This  result  for  y  (x,p)  follows  by  induction  on  n„ 

n 

Finally,  the  uniform  convergence  of  the  series  (A. 11)  gives  the  result 
that  y(x,p)  is  uniformly  continuous  in  R. 

Theorem  2;  If  H(G,z)  Is  an  analytic  function  of  its  complex  arguments 

in  the  region  R  of  G  — -  z  space  given  by 

|z  -  Zq|  <  h  ,  |G  -  G0|  <  Mh  , 

and  if  conditions  |H(G,z)j  <  M  and  |H(G,z)  -  H(G*,z)|  <  NjG~G*|  are 

imposed  on  H  in  R,  then  the  solution  G(z)  of  the  differential  equa- 
cLG 

tion  —  =  H(G,z)  which  satisfies  the  condition  G(zq)  =  Gq  exists  and 
is  a  unique  analytic  function  of  the  complex  variable  z  in  the  region 

R. 


-  85  - 


Proof;  Since  R  is  an  open  region,  here  no  question  of  uniform  contin¬ 
uity  arises.  The  proofs  of  existence  and  uniqueness  go  through 
as  in  Theorem  1  with  x  and  y  replaced  by  z  and  G  respectively  and 
with  the  appropriate  replacement  of  <  with  <  due  to  the  different  def¬ 
inition  of  R  .  It  is  assumed  that  paths  of  integration,  where  they  occur 
lie  inside  R. 

We  now  prove  the  analyticity  of  G(z).  The  equations  correspond¬ 
ing  to  (A. 07)  and  (A.08)  are 

°o(z>  ■  Go 

and  Gn(z)  *  Go  +  f  HtGn  i(2)»z^dz  • 

Jzo 

Since  the  integral  of  an  analytic  function  is  analytic,  then  G^(z)  and, 
by  induction  on  n,  ®n(z)  are  analytic  in  R.  Again  the  uniform  con¬ 
vergence  of  the  series  (A. 11)  ensures  the  analyticity  of  G(z). 

These  results  are  valid  when  the  boundary  condition  is  satisfied 
at  a  point  y^  inside  the  region  R.  They  also  hold  if  the  condition  is 
met  at  a  point  on  the  boundary  as  is  shown  in  Graves , "Theory  of  Functions 
of  Real  Variables,  Second  Edition",  McGraw-Hill,  1956,  page  157* 
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